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ABSTRACT 

This thesis investigates various techniques for fore- 
casting the volume of containerizable cargo that flows into 
the container stuffing station at the Military Ocean 
Terminal, Bay Area, Oakland, California. Cargo input data 
is analyzed in terms of weekly cargo volume inputs for a 
selected number of major ports of debarkation. The time- 
series data for these ports is first tested for serial 
correlation. Based on the affirmative results of the 
serial correlation test, the following forecasting methods 
are investigated: the moving average, the exponentially 

weighted average, the exponentially weighted average \^ith 
trend adjustment and the exponentially weighted average with 
an adaptive response rate. By means of statistical testing 
procedures, the "best" forecasting method is determined. 
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I. INTRODUCTION 



In light of the dispersal of United States forces 
throughout the world, the shipment of the vast quantities 
of needed supplies and material to overseas stations is a 
predominant logistics problem for the Department of Defense. 
Since the late 1960's the need for efficient and economical 
handling of this cargo has led the DoD to use containerized 
shipments on commercial sea-going carriers as the primary 
routine transportation source. The retention of break bulk 
type shipping operations has been necessary due to the size 
limitations of containers. Also, the fulfillment of urgent 
requirements where time is critical have been met using air 
transport. However, ocean containerization offers the most 
cost-effective method of handling routine shipments between 
ports equipped to handle this type of traffic, as are most 
Pacific Ocean ports. 

The ability to fill or stuff a container at the con- • 
signor's warehouse and unpack it at the consignee's ware- 
house, with no intermediate handling of the cargo, is a 
tremendous asset. However, many shipments are not large 
enough to economically warrant their own vans; therefore, 
stuffing with other shipments bound for the same destination 
may be accomplished at a central location. One such stuff- 
ing facility is the container freight division (CFD) of the 
Military Ocean Terminal Bay Area (MOTBA) located in Oakland, 
California. This facility stuffs vans with containerizable 
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cargo enroute to ports, mainly in the Pacific and the Far 
East, and arranges transportation of the containers to the 
ports via commercial carrier. 

All commercial carriers require that bookings for vans 
and space aboard vessels be made in advance of the ship's 
loading date. As of the fall of 1974, this lead time was 
at least three weeks. To effectively make these bookings 
for an advanced date, a prediction of the cargo on hand on 
the date of arrival of the vans is required. Normally, 
vans arrive at the stuffing station four to six days prior 
to the loading date. Thus, the request for space must take 
into account Cl) the present volume of cargo in the ware- 
house, (2) the forecasted arrivals of cargo during the 
period between the booking and the van arrivals for stuff- 
ing and (3) any previous bookings that will deplete the 
cargo volume destined for the port during the forecasting 
period . 

The general objective of this thesis is to examine 
practical forecasting methods which can be utilized to pre- 
dict cargo arrivals at the container stuffing station. To 
predict arrivals there are two general approaches that 
could be used. The first approach is based on notification 
of shipment departures by consignors, which would provide 
advanced notice of arrivals at the stuffing station. The 
other approach is based on historical data, from which pro- 
jections for future arrivals may be drawn. 

An advance notification system does not presently exist 
for the majority of cargo arriving at the station; however, 
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it could be produced by requiring consignors to provide the 
CFD with three-week advance notification of cargo shipments 
perhaps by post card. Use of this approach would require 
analysis of the delay between the arrival of notifications 
and shipments from each consignor. Also, the statistical 
characteristics and stability of the transit times for ship 
ments and notifications must be determined. It is doubtful 
that most consignors have enough advance knowledge of a 
shipment's departure to provide a timely and accurate noti- 
fication, especially if the departure is two or three weeks 
in the future. The effect of the possible lack of consigno 
cooperation must also be assessed. Taking into account the 
resources that would be required to effectively implement 
an advance notification approach, it was decided that fore- 
casting methods which rely only on historical data for 
future projections would be more appropriate for a thesis 
topic. 

The historical data used here consisted of the cargo 
arrival data for the CFD for fiscal year 1974; it repre- 
sented the most recent peacetime data available. The data 
provided the volume, weight and date of arrival of each 
shipment. Of the three alternative units for forecasting 
(number of shipments, volume and weight), the volume of the 
cargo was selected because it is currently used by the mari 
time community as the measure for booking containers on 
vessels . 

Two alternative aggregation levels could have been used 
for forecasting. One was based on the port of debarkation 
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(POD), and the other based on consignee destination of the 
cargo. There are more than 1,000 active consignees served 
by the CFD, which would require a very large number of in- 
dividual forecasts if the consignee base were used. PODs 
are the currently used aggregation levels for forecasts by 
CFD personnel because van space is, by convention, booked 
aboard ships destined for PODs, not consignees. Therefore, 
PODs were selected as the aggregation level. Analysis of 
the data revealed 16 major PODs which received over 90% of 
the shipments passing through the CFD. 

The effectiveness criteria used for the selection of an 
optimal method are: 1) computational feasibility, and 2) 

the ability of the method to maximize forecasting accuracy 
by minimizing forecasting error. Computational feasibility 
is limited here to those methods which do not require the 
use of a computer to routinely generate forecasts. Mean 
absolute deviation and mean deviation squared are used as 
measures of forecast errors. 

Forecast errors were computed over time horizons of 
one, two and three weeks, corresponding to possible lead 
times in the booking process. In other words, for a two- 
week time horizon a single forecast was used for two suc- 
cessive weeks. The deviation for the forecast was computed 
by taking the average difference between the forecast and 
the actual input volume for each of the weeks covered by 
the forecast. 

Generally, the selection among various methods of fore- 
casting depends upon the statistical characteristics of the 
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data. These characteristics are differentiated by the degree 
of statistical dependence between successive observations. 

If there is no statistical dependence between successive 
observations, then ordinary trending and indexing techniques 
may be applied. On the other hand, the presence of statis- 
tical dependence among successive observations, (called 
serial correlation) , implies a need for methods which place 
more emphasis on recent past observations when predicting 
future volumes. 

Tests for serial correlation of the data for the major 
PODs were made using the Durbin-Watson d-statistic. The 
aggregated results of these tests indicated that the pre- 
sence of serial correlation was statistically significant. 

Figure 1 is a graph of the weekly input volumes of POD 
UL7 for a 20-week period commencing with the nineteenth 
week of FY '74. This was selected as a typical representa- 
tion of the data for the major PODs. This particular POD 
exhibited strong positive serial correlation. 

With the existence of serial correlation in the data 
confirmed, the moving average technique and three exponen- 
tially weighted average techniques were examined for pre- 
diction accuracy. These methods were chosen for their 
relative ease of computation and ability to adapt to changes 
in volume inputs over time. 

Simple moving averages [Ref. 5] of from one to ten week 
durations were the first methods explored. In each method, 
the most recent weeks for the designated duration are aver- 
aged and the result applied to the forecast through the 
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Figure 1: Weekly volume inputs for POD UL7 



designated horizon. The three-week moving average was the 
technique in actual use at the CFD. The resulting data in- 
dicated that for all time horizons, moving averages greater 
than three weeks generally provided lower deviations and 
thus, better forecasts than the three-week moving average. 
Note that the longer the duration of the average, the less 
responsive the predictions are to recent changes. The data 
indicated that over-sensitivity was as undesirable as no 
sensitivity in most cases. 

The next method examined was the simple exponentially 
weighted average, [Ref. 5] , or exponential smoothing method. 
This method used a weighting or smoothing constant, alpha, 
a value usually between zero and one, to regulate the em- 
phasis on the most recent observation. Large values of alpha 
place increasing emphasis on the most recent observations 
and less on preceding observations. A forecasted average 
for an observation period was computed as the sum of the 
observations in that period times alpha, plus the forecasted 
average for the previous period times one minus alpha. In 
this way, the new observation was added to the accumulated 
weighted average, a continually changing number. This fore- 
cast average was then used as the forecast and applied over 
the next weeks of the designated horizon and the deviation 
statistics accumulated. 

To discover the levels of alpha which minimized error 
measurement, the levels of alphas were varied by 0.01 from 
-0.1 to 1.0, and applied to the data for each POD and hori- 
zon. The alpha levels which produced minimal deviation 
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statistics ranged from -0.02 to +0.77 depending on the POD 
and time horizon specified. The optimal single alpha level 
for all PODs was about 0.35 and varied slightly depending 
on the time horizon. 

The next forecasting method explored was the exponen- 
tially weighted average with adjustment for trends [Ref. 5], 
This procedure was examined to determine if trend adjust- 
ments would reduce forecast error. This procedure started 
with the same basic formula as the simple exponentially 
weighted average, however a trend factor is also computed. 
This trend is accumulated and smoothed in the same manner 
as was applied to the forecast averages, resulting in an 
average trend adjustment. The forecast volume to be used 
in the predictions was created by combining the forecasted 
average with the average trend adjustment. 

To obtain minimal errors, alpha levels were again varied 
by 0.01 from -0.10 to 1.00 and the method applied to the 
data for all PODs. The results revealed minimized devia- 
tions for alpha levels between -0.01 and 0.38, again, de- 
pending on the POD and time horizon specified. When 
accounting for trends, the optimal single alpha level for 
all PODs was 0.16, varying slightly for the different time 
horizons . 

The final forecasting procedure used was the exponen- 
tially weighted average using an adaptive response rate 
technique to modify the alpha level used from period to 
period to account for fluctuations in the observations 
[Ref. 18]. The absolute and actual errors, or deviations 
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and actual errors between the observed and forecasted 
volumes were separately accumulated and smoothed in the 
same fashion as the forecast average using a smoothing con- 
stant, gamma. A tracking signal was then computed by di- 
viding the smoothed absolute error by the smoothed error. 
Deviation statistics using the simple exponentially weighted 
average were computed with the alpha level for each forecast 
set at the absolute value of the tracking signal. 

To evaluate the adaptive response rate method, it was 
applied to the data for each of the PODs while varying the 
time horizon and the level of the smoothing constant, gamma, 
by 0.01 from -1.00 to 1.00. Here the use of negative con- 
stants indicates that the new smoothed error term is created 
by enlarging the old smoothed error and adding or subtract- 
ing a portion of the new error to it. The minimum deviation 
statistics for the individual POD and various time horizons 
were realized with gamma levels ranging from -1.00 to 0.49. 
The resultant optimal single gamma level for all POD was 
about -0.05 depending on the horizon used. 

In summary, the experimentation with the above forecast- 
ing methods had created deviation statistics for the follow- 
ing forecasting procedures: 

1) Three-week moving averages (maintenance of status 
quo) ; 

2) Moving averages with a single, aggregately optimal 
duration ; 

3) Moving averages with the optimal duration for each 

POD; 



20 



4) Simple exponentially weighted averages using a 
single, aggregately optimal alpha level; 

5) Simple exponentially weighted averages using in- 
dividually optimal alpha levels for each POD; 

6) Exponentially weighted averages with trend adjust- 
ment using a single, aggregately optimal alpha level; 

1 ) Exponentially weighted averages with trend adjust- 
ment using individually optimal alpha levels for each POD; 

8) Exponentially weighted averages with adaptive 
response rates using a single, aggregately optimal gamma 
level ; 

9) Exponentially weighted averages with adaptive re- 
sponse rates using individually optimal gamma levels for 
each POD. 

Using these deviation statistics, the selection of the "best" 
forecasting method was made subject to the following 
considerations: 

1) For operational consistency, it was assumed that 
only one of the methods would be chosen for use with all 
PODs . No attempt was made to evaluate the problem by 
choosing the optimal forecasting method for each POD. 

2) The deviation statistics for the one week time 
horizon were used in the evaluation as the relative effec- 
tiveness of the various methods for each POD. 

3) The mean absolute deviation statistic was selected 
as the principal goodness criterion since no extra penal- 
ties for large individual forecast errors were intended. 
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In the end, statistical significance tests are used to 
reveal whether the mean absolute deviations using one fore- 
casting method were significantly less than the deviations 
of another method. The tests were conducted at a signifi- 
cance level of 0.05. The results indicated that the simple 
exponentially weighted average and adaptive response rate 
methods were not significantly different from each other. 
However, both methods were better than any of the others. 

Based on the above results the use of the simple expo- 
nentially weighted average was chosen. The selection was 
founded on the ease of computation of the forecast, as the 
adaptive response rate method required many more steps for 
each calculation and offered no better results. Addition- 
ally, the adaptive response rate method required a broad 
range of gamma levels and was characterized by multi-modal 
relationships between the gamma level and MAD for most cases 
of PODs and time horizons. This indicated great sensitivity 
of accuracy to the gamma levels for particular data sets. 

When compared to the three-week moving average, the use 
of this simple exponentially weighted average technique 
provided 4.9%, 8.0% and 10.0% reductions in mean absolute 
deviations for the one-, two- and three-week time horizons 
respectively. In addition to increasing the forecast ac- 
curacy the exponentially weighted average is also as easy 
computationally as the three-week moving average. Opera- 
tional benefits of the increased forecast accuracy are a 
reduction in the average age of cargo and volume of cargo in 
the warehouse, and increased cube utilization. 
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II. BACKGROUND 



A. THE BEGINNING OF CONTAINERIZATION 

While containers had been used for transporting cargo 
aboard vessels as early as the 1820's, the modern era of 
containerization began in 1966 when Malcolm McLean of Sea- 
Land Service, Inc., began loading containers aboard ships 
bound for Europe. The economics of containerization soon 
became apparent. By the end of 1968, one-fourth to one- 
third of shipments on major U.S. routes were containerized. 
American-f lag operators as well as their foreign-flag com- 
petitors began to invest heavily in the containerized sys- 
tem. This system includes not only containers, but also 
specialized ships, container terminals and road equipment. 

Within the industry the container came to be called the 
greatest thing in packaging since the paper bag. Most con- 
tainers are simple boxes of plywood or aluminum and are 
manufactured in standard sizes: approximately 8 feet wide, 

8 feet deep and various lengths between 20 to 40 feet long. 
The container volume capacity ranges from 1,100 cubic feet 
to 2,400 cubic feet, and weight capacity ranges from 38,000 
pounds to 46,000 pounds depending on the size of the van. 

Containerization required the shipping industry to 
adopt an entirely new line of thinking and provided the in- 
dustry with many phases of economic improvement: more ef- 

ficiency in cargo and terminal handling, a reduction in 
merchandise damage due to weather and handling, a reduction 
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in terminal facility requirements, a savings in reduced 
special packaging for export, a reduction in pilferage and 
other losses, the ability of manufacturers to operate with 
reduced inventory, and the benefit of computerized inven- 
tory-control systems. 

Containerization has cut shiploading and unloading time 
from as much as eight days to as little as twelve hours. 

The varied forms of containerization are enabling cargo 
ships to spend as many as nine out of every ten days' pro- 
ductivity at sea instead of half their time in port -- 
causing some shipping lines to have alternating crews so 
that the men can spend some time with their families. A 
typical containership is loaded with some 700 vans, seven 
deep and six across, and is on its way at 22 knots. 

The new era has resulted in transmodal systems to move 
cargo via land, sea and air with computerized efficiency. 

In some cases, overall transport costs are down 30-40% as 
compared to conventional break bulk cargo handling. 

Sea-Land Service, Inc., has become the biggest container 
operator in the world and is one of the two largest U.S. 
steamship lines. Sea-Land has five 33-knot containerships 
with capacity for 1,082 units. With other carriers taking 
advantage of container economics, the United States flag 
carriers have become the world's largest containership 
fleet . 

The port of Oakland, California is now the leading West 
Coast cargo shipping terminal -- mainly due to containeriza- 
tion. The port handles 45% of all container shipping on 
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the West Coast and is the second largest container port in 
the world -- New York being the largest. 

According to James L. Reynolds, President of the 
American Institute of Merchant Shipping (AIMS) , by 1976 
cargo moving to and from U.S. ports is expected to rise to 
270 million measurement tons (a unit of volume representing 
40 cubic feet) plus 7.5 million in military cargo and an- 
other 20.5 million in government impelled cargo such as 
foreign aid shipments. 

Containerization has been recognized as a system that 
will best serve the needs of the public and the transporta- 
tion industry. The container revolution stems from the im- 
mense transportation capability gained from moving goods 
swiftly, safely and reliably. 

B. CONTAINERIZATION AND THE DEPARTMENT OF DEFENSE 

In recent years the United States Department of Defense 
has used commercial containerization as the primary means 
of transporting its general ocean cargo. Substantial cost 
savings from the use of containerships in lieu of using the 
traditional break bulk cargo ships have been recognized. 

Yet, there still remain problems for transportation offi- 
cials within DoD to reckon with. 

Figure 2 illustrates the portion of the transportation 
system within the Department of Defense to which this thesis 
addresses itself. 

Orders are placed by DoD overseas commands known as 
consignees. Various vendors and military depots known as 
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Figure 2: Export Cargo Movement 




consignors, fill these purchase orders which the consignees 
need to execute their missions. It now becomes the job of 
the DoD transportation managers to deliver these items to 
the consignee. 

Cargo leaving the consignor for the consignee may be 
transported across the ocean in either container or break 
bulk ships. Conventional ships transport break bulk cargo 
to the requesting overseas command. If the goods are 
ordered in large enough quantities and are containerizable , 
the consignor may place the cargo in containers at his ware- 
house. This is known as source stuffing. The sealed con- 
tainer is then shipped to the POE where it has been booked 
aboard a containership for overseas transport to the con- 
signee. This type of operation provides the most rapid and 
economical form of ocean shipping. 

Another situation occurs when container stuffing sta- 
tions, at the POE or inland, receive quantities of contain- 
erizable cargo for consolidated stuffing into containers. 

At the stuffing station, cargo for particular consignees is 
stuffed only after meeting volume and weight restrictions. 

A commercial shipping company then receives the containers 
from the stuffing station for lift aboard containerships 
and overseas transport. 

At the Port of Debrakation (POD) single consignee con- 
tainers are sent directly to the respective consignees. 
Containers with cargo for more than one consignee are trans- 
ported to designated break bulk stations (BBS) , which may 

itself by a consignee, where they are broken down for furthe 
shipment to individual consignees. 
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C. WESTERN AREA, MILITARY TRAFFIC MANAGEMENT COMMAND 



Western Area, Military Traffic Management Command (WAMTMC) 
is a jointly-staffed field organization under Headquarters, 
Military Traffic Management (MTMC) , Washington, D. C. WAMTMC 
is responsible for transporting export cargo of the Depart- 
ment of Defense and other government agencies to overseas 
commands. It has been located at the Oakland Army Base, 
Oakland, California, since being established in February 
1965. MTMC is a major field command of the Department of 
the Army. 

WAMTMC is responsible for transportation management of 
domestic and export shipments in 14 western states: Idaho, 

Arizona, California, Montana, Nevada, Oregon, Utah, New 
Mexico, Washington, Colorado, Nebraska, South Dakota, North 
Dakota and Wyoming. One of the three Military Ocean Termi- 
nals operated by WAMTMC on the West Coast will receive ex- 
port cargo from these 14 states which is destined for 
shipment to overseas installations. The three terminals 
operated by WAMTMC are: Military Terminal Unit, Pacific 

Northwest, Seattle, Washington (PNW) ; Military Ocean Termi- 
nal, Bay Area, Oakland, California (MOTBA) ; Southern 
California Outport, Long Beach, California (SCO). During 
peacetime, WAMTMC also directs the flow of all air cargo 
in the Military Airlift Command's transportation system 
through the Military Airlift Clearance Authority (MACA) . 

D. MILITARY OCEAN TERMINAL, BAY AREA 

Military Ocean Terminal, Bay Area is the consolidation 
of the Army and Navy terminal facilities in the San Francisco 
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Bay Area. It represents the largest operating unit of 
WAMTMC with terminal facilities at the Oakland Army Base 
and the Alameda Reefer Facility. Seven deep water berths 
are under its control as well as an 84-acre tidewater con- 
tainer stuffing area located at the Oakland Army Base. 

Also under MOTBA's control are the ports of Stockton, 
Sacramento, and Eureka. 

E. TIDEWATER CONTAINER STUFFING STATION 

The Container Stuffing Station (CSS) operates under the 
Container Freight Division (CFD) at MOTBA. Containerizable 
cargo is loaded into the shipping containers at both instal- 
lations. These containers are then transferred to the com- 
mercial shipping companies for shipment to overseas commands. 

The CSS is operated under civilian contract. Container- 
izable cargo is placed in warehouses until containers become 
available for stuffing. In most cases commercial containers 
are used for stuffing, although military vans are sometimes 
used. The contractor has certain restrictions on container 
stuffing. Some examples are cargo mix, container utiliza- 
tion, and consignee mixing limitations. In addition, volume, 
weight and size of containers vary among the various com- 
mercial carriers. There are 16 types of containers presently 
in use; and, in general, each company can accommodate only 
those containers which have been specifically designed for 
its vessels. 
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F. MILITARY SEALIFT COMMAND, PACIFIC 



The Military Sealift Command *(MSC) is an agency that 
coordinates the various ocean transportation activities of 
the Department of Defense. Within the continental United 
States (CONUS) the relationship between the Military Sealift 
Command, Pacific (MSCPAC) and WAMTMC is very close. In 
order to get its cargo aboard break bulk or container ves- 
sels, WAMTMC notifies MSCPAC of the amount of space required, 
the number of containers needed by the shipper and the size 
of the containers. MSCPAC books the space aboard a vessel 
and notifies WAMTMC of the vessel sailing date, lift date 
for the containers and the number of containers booked. 

This procedure is carried out whether the shipper is a ven- 
dor or WAMTMC is the shipper. 

G. MTMC MANAGEMENT INFORMATION SYSTEMS 

MTMC employs two management information systems in order 
to control the flow of cargo within CONUS that will be 
eventually shipped to overseas commands. Both systems 
utilize a Burroughs 5500 computer. The two systems are the 
Mechanized Export Traffic System (METS) and the Surface Ex- 
port Cargo System (SURS) . 

METS monitors the DoD cargo flow from its origination 
at the consignor's warehouse until its arrival at a state- 
side Port of Embarkation, i.e., one of WAMTMC' s three ocean 
terminals. The consignor, which may be a commercial vendor 
or a military depot, thus provides the basic input into the 
system. This input includes shipment data, commodity 
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information, transportation data, destination, etc. Im- 
properly formated documentation is frequently received from 
the commercial vendors requiring the correction of these 
errors by manual methods. 

Once cargo is within one of the three west coast ocean 
terminals, it falls under SURS, or SURS/CARDPAC . A CARDPAC 
is a set of eight IBM cards containing various shipping 
data. These cards have been punched by terminal personnel 
with data obtained from METS and any additional information 
that may be needed from shipping forms. 

The CARDPAC; the Expected Receipt List (ERL), a list 
of the cargo expected in a particular shipment; and the 
applicable source document, e.g., a Government Bill of 
Lading (GBL) , Commercial Bill of Lading (CBL) , Dray Tag, 
or Transportation Control and Movement Document (TCMD) are 
forwarded to the terminal's cargo receiving area. As the 
cargo is unloaded a check is made against information con- 
tained on the CARDPAC and corrections are made if necessary. 
The CARDPAC is then attached to the shipment and will be 
utilized to monitor its flow through the terminal facilities. 

In summary, by using the METS and SURS management infor- 
mation systems, WAMTMC is able to monitor the flow of cargo 
into and within its terminal facilities. Management reports 
generated by the system and statistical analyses based on 
the information provided by the system are only as good as 
the data going into it. 
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H. CARGO CLASSIFICATION 



There are two types of cargo passing through the Con- 
tainer Stuffing Station: (1) release unit (RU) ; and (2) 

less than release unit (LRU) . 

Shipments in excess of 10,000 pounds and special cate- 
gory cargo, such as classified cargo, are classified as 
release unit material. Under DoD Regulation 4500. 32-R, 
Military Transportation and Movement Procedures (MILSTAMP) , 
release unit material requires positive export traffic 
release . 

Consignors having break bulk or containerizable cargo 
that fall into RU classification must request and receive 
clearnace to move that material before it can be shipped. 
WAMTMC , Export Control Division, serves as the clearance 
authority for RU material for the fourteen Western states. 

This thesis is not concerned with break bulk or source 
stuffed cargo. Therefore, this paper will not go into a 
detailed description of the booking procedure except to 
mention that it is accomplished through the coordinated ef- 
forts of MSCPAC, WAMTMC, the shipper and the commercial 
carrier . 

It should be mentioned that if RU classified cargo is 
not of sufficient quantity to be source stuffed, it is 
shipped to the CFD for stuffing. 

Eighty percent of the shipments passing through MOTBA 
are of the LRU classification. LRU material does not re- 
quire positive release into the system. Therefore, there 
is an uncontrolled flow of cargo coming into the Container 



32 



Freight Division. This results in a very difficult booking 
problem for the Container Freight Division, MOTBA, WAMTMC 
and MSCPAC. It is the forecast of this uncontrolled cargo 
flow to which this thesis addresses itself. 

Figure 3 represents the cargo cycle to be considered. 

The procedures at each step are as follows [Ref. 15]: 

1) A DoD overseas command makes a request for goods 

by requisition. The specific requestor is the user of these 
goods and the ultimate consignee. The requisition is to be 
filled by the shipper (consignor) who may be a commercial 
vendor or a military depot. 

2) Material that is classified as less than release 
unit is shipped by the consignor directly to the Container 
Freight Division at MOTBA. 

3) Under Military Traffic Management Regulation, con- 
tainer requirements are offered to WAMTMC by MOTBA. 

4) A request for booking is then submitted to MSCPAC. 

5) MSCPAC then books space aboard ocean carrier. 

6) Once booking aboard the vessel is made, the related 
information is forwarded to WAMTMC. 

7) The Container Freight Division at MOTBA receives 
the booking data from WAMTMC. 

8) The Container Freight Division communicates directly 
with the ocean carrier to arrange for spotting and pick-up 
of containers. 

9) Once containers have been stuffed at the CSS, they 
are delivered directly to the ocean carrier's container yard. 
From this point shipment from POE to POD is made. 
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Figure 3: Supply and Transportation Cycle for Containers [Ref. 15]. 




10) Delivery of material is made to the consignee by 
either the ocean carrier or by arrangements made by the 
overseas command. 

I. CURRENT BOOKING PROCEDURES 

As previously stated, a request for booking is made 
twenty-one days in advance and is based on a forecast of 
the volume of cargo in measurement tons that is expected to 
be on-hand at the time the vans are to be stuffed. An of- 
fering is made based on the following formula: 

Booking _ ( Cargo- ^ _ , Cargo + , Forecasted ^ 

Offering ^ on-hand * ^ Booked ' Cargo Receipts * * 

The CSS computes the previous three weeks average weekly 
receipts and uses this figure as a basis for forecasting the 
volume of cargo arriving at the facility between the date of 
offering and the stuffing of containers. Figure 3 illustrates 
the present method. 

Figure 4' depicts the situation where a particular POD 
has 100 measurement tons (MT's) of cargo on the floor on 
Julian calendar day 70. To keep the computations simple, 

250 MT's of cargo will represent the average weekly receipts 
for this POD for the past three weeks (15 days, since only 
weekdays are included) giving average daily receipts of 50 
MT's. Container space has already been booked aboard three 
vessels during the 21- day period under consideration. 

The cutoff dates have been made for days 77, 79 and 84 for 

360, 80 and 120 MT's respectively. Using the past three 

♦ 

weeks average receipts as the basis of the forecast, there 
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in MTs 400 
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Figure 4: Booking Procedure 



will be 360 MT’s on-hand on day 77. The two-day period 
where no cargo accumulation occurs represents a Saturday 
and a Sunday. The five working days will accumulate 250 
MT's. Since space for 360 MT’s has been previously booked, 
all cargo on the floor for this POD will be stuffed. Two 
days later there will be 100 MT's on the floor but 20 MT’s 
will not be stuffed. On day 84 only 120 of the 170 MT’s 
will be stuffed. By day 91, the day for which the forecast 
is made, the forecast says there will be 300 MT’s on the 
floor for this particular POD. The number of vans requested 
will depend on the average cube utilization for the POD. 

If the POD's average is 50 MT's then an offering for six 
vans will be made on day 70 to be available for stuffing 
on day 91. 

J. CANCELLATION PROCEDURE 

Cancellation of space booked aboard a carrier must be 
made no later than "a reasonable length of time" prior to 
the cutoff date according to MSC’s Container Agreement and 
Rate Guide, RG8. If cancellation cannot be made within this 
time and if the carrier cannot utilize this space, then the 
government may be charged for the space even though it was 
not used. The present policy at WAMTMC is to cancel no 
later than five days prior to cutoff. 

K. OPERATIONAL RESTRICTIONS AND REQUIREMENTS 

It should be noted that there are other operational re- 
strictions and requirements that affect stuffing and booking 
procedures, although a thorough understanding of them is not 



essential to this work. Briefly they are: 1) a first-in- 

first-out (FIFO) cargo stuffing procedure requiring a mini- 
mum average FIFO performance of 80% for each POD, 2) shipment 
priorities that are exceptions to the FIFO procedure, 3) a 
requirement for the timely handling of certain goods such 
as household effects, 4) a requirement for a monthly cube 
utilization average of 75% with no less than 50% cube utili- 
zation to ensure that container utilization is more economi- 
cal than break bulk shipping, 5) a restriction on cargo 
mixing such as Military Assistance Program (MAP) cargo not 
being permitted to be mixed with any other type of cargo, 

6) the requirement of contracting the low cost carrier if 
there is more than one that can provide the required ser- 
vice, and 7) a ceiling of 4,000 measurement tons on the floor 
at any one time. 



38 



III. STUDY OBJECTIVES AND PROCEDURES 



A. STUDY OBJECTIVES 

The essence of this study is to observe the flow of cargo 
into the Container Freight Division at MOTBA and attempt to 
forecast by statistical and analytical techniques the expected 
cargo inputs at some future point in time. The forecast is 
of the volume of containerizable cargo passing through the 
CFD destined for a particular POD. Break bulk cargo volume 
and source stuffed vans are not considered in this analysis. 

The necessity for such a forecast evolves from the fact 
that the Military Sealift Command must book container space 
aboard vessels as much as three weeks in advance of lift. 
Because of the demand for container shipping in recent years, 
commercial carriers are able to make this demand on their 
customers. It remains to be seen whether a shift in the sup- 
ply and demand factors will alter this requirement to any 
significant extent. 

Although forecasting and forecasting models are generally 
associated with inventory and production control, the same 
general principles apply to the situation at MOTBA. Es- 
sentially, the situation is one of demand -- the demand by 
overseas commands for material support. In this respect the 
forecasting problem can then be broken down into components 
of that demand. Examples of these components are average 
demand, trend effects, seasonal effects, and noise or random 
effects. By recognizing these components one can then con- 
struct statistical models to deal with them. 
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Accurate forecasting may have an influence on some of 
the performance variables that the container booking system 
generates. These are frequency of cancellations due to 
overbooking, age of cargo at stuff, cube utilization, and 
total volume of cargo on the floor. 

In forecasting, the past is projected into the future. 

A good forecasting model should provide information feedback 
accurately and quickly so that the manager can make the 
necessary and appropriate adjustments. It should provide a 
means of reflecting sudden shifts caused by policy changes 
within the system or uncontrollable exogenous factors which 
may also influence the system. A forecast that is based on 
out-dated information will not be beneficial to a manager 
who must book container space three weeks in advance. 

A forecast is only as good as the data it uses. In this 
thesis it is assumed that the historical cargo flow data 
used to evaluate the forecasting methods is accurate. Data 
totals obtained from the computer tapes were verified against 
the totals that were manually computed at MOTBA. 

B. PROCEDURES 

As stated previously the system variable that this thesis 
addresses is the volume of containeri zable cargo for a POD 
flowing into the Container Freight Division at MOTBA. Al- 
though there are weight restrictions to consider, these re- 
strictions have not been significant factors in the booking 
procedure . 
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The data used throughout this thesis was taken from 
actual operations at the Container Stuffing Station. The 
data covers an 18-month period dating from 1 January 1973 
to 30 June 1974. This period generated an enormous volume 
of data which reflected the magnitude of the container 
stuffing operations at MOTBA. 

In order to handle this volume of information the ana- 
lysis was done in four separate phases. These were: 

1) a distribution audit in order to identify the vari- 
ables needed in the analysis, 

2) the grouping of data in a more usable format, 

3) a test for serial correlation in the data, and 

4) the application of various forecasting techniques 
and their evaluation. 

In performing the distribution audit some of the com- 
puter programs from the McCarthy and Carter thesis were 
utilized. Only slight modifications to these programs were 
needed to manipulate the enormous volume of data generated 
by the eighteen months of operations. The audit identified 
the number of shipments and the volume of these shipments. 
From this information one could identify major Ports of 
Debarkation, thus narrowing the scope of the problem. 

In order to arrange the data into a more usable format 
a number of short computer programs were written. These 
programs grouped the data from daily to weekly totals in 
the aggregate as well as POD levels. 

Since this work dealt with time-series analysis it was 
reasonable to test the data for serial correlation in order 
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to limit the range of potential forecasting methods. A 
computer program was utilized to calculate the statistical 
information needed to evaluate the various PODs. 

Several computer programs were then used to apply various 
forecasting techniques to the data. The simplest techniques 
were applied first. In some cases these initial techniques 
were expanded to take into consideration the components of 
demand previously mentioned. The various forecasts were 
then evaluated based on the difference between the fore- 
casted and actual volume. 

All the computer programs used in the data analysis were 
written in the FORTRAN IV computer language. These programs 
were run on the IBM 360/67 computer at the W. R. Church 
Computer Center, Naval Postgraduate School, Monterey, 
California . 
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IV. DATA ANALYSIS 



A. DISTRIBUTION AUDIT 

Raw data was taken from a set of 19 computer tapes sup- 
plied by WAMTMC and was combined with existing data from 
previous research work to form a set of data covering 18 
months. The period covered was Quarters III and IV of FY 
*73 and all four quarters of FY '74. Computer programs 
were used to extract the pertinent information from the 
tapes, combine it with existing data and place the aggre- 
gate information in alphabetical order by POD and consignee 
in a data file. 

Computer Program Two of the McCarthy and Carter thesis 
was used to delineate the data for further analysis. This 
program determined the number of shipments, the total volume, 
the average volume, the standard deviation of the volume, 
the total weight, the standard deviation of the weight, the 
average density and the standard deviation of the density 
for each POD. The program also grouped each POD's shipments 
into intervals of volume, weight and density. 

The next procedure entailed the examination of the gross 
weekly totals of shipments and volume of cargo arriving at 
the container stuffing station during the 18-month period. 

The following information was obtained: 



SHIPMENTS 



Quarters III 8 IV, FY '73 
FY '74 

18-month Totals 



111,463 shipments 
124,592 shipments 
236,055 shipments 
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VOLUME 



Quarters III $ IV, FY '73 
FY '74 

18-month Totals 



8,398,931 cubic feet 
10,070,724 cubic feet 
18,469,655 cubic feet 



There was obviously a tremendous amount of activity at 
MOTBA during the second half of FY '73. Over 47% of the 
total shipments and 45% of the total volume for the 18-month 
period occurred during the first six months of the data set. 

It was concluded that residual Vietnam conflict effects 
were responsible for generating this high volume of activity 
during the first six months of the 18-month period and that 
it did not represent the present and the expected future 
flow of cargo at MOTBA. Therefore, it was decided to use 
only the data for FY '74. 

The audit was run a second time using only FY '74 data. 
During that period 157 PODs and 3,356 consignees were 
identified as having received shipments during this period. 
These PODs and consignees accounted for 124,592 shipments 
with a total volume of 10,070,724 cubic feet as previously 
stated. This was an average weekly total of 2,396 shipments 
with an average weekly volume of 193,668 cubic feet (or 
4,842 measurement tons). Total weight for the period was 
176,213,333 pounds. 

In reviewing the yearly gross totals for the 157 PODs, 
it was obvious that 10% of these dominated the activity at 
MOTBA. Sixteen PODs were identified as being the most ac- 
tive in terms of volume of cargo arriving at the container 
stuffing station during the period. The cutoff point was 
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arbitrarily chosen at 100,000 cubic feet of volume per year. 

As a result the least active of the 16 major PODs averaged 
a volume input of 2,041.35 cubic feet per week (or approxi- 
mately 51 measurement tons). The number of shipments for 
these major ports totalled 115,732 or 92.89% of all ship- 
ments for the period. Their volume total of 9,370,693 cubic 
feet was 93.05% of the total volume. 

Table I, Appendix A, shows the weekly gross totals of 
shipments, volume and weight arriving at the CSS in the second 
half of FY '73. The totals for FY '74 appear in Table II. 

Table III, Appendix A, shows the data summary of the dis- 
tribution audit. Tabular data for the 16 major PODs is 
presented in Table IV. Table V lists the activities of the 
major PODs by shipments, volume and weight. 

The computer output of the intervals of weight into 
which each of the shipments fell provided a look at the por- 
tion of the release unit shipments arriving into the system. 

If the number of release units were significantly large, it 
may be possible to forecast with greater accuracy the volume 
of cargo expected to arrive during a specified period. 

There were 1,900 shipments of ten thousand pounds or 
more during FY '74 using only data generated by the 16 major 
PODs. These shipments should have been classified as re- 
lease unit shipments and it is assumed that they were. They 
represented 1.64% of the shipment totals for the major ports. 

It was concluded that the amount of release unit ship- 
ments arriving at the container stuffing station was insignifi- 
cant in view of these figures. Furthermore, operations at 
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MOTBA indicate that the expected arrival date of RU cargo 
was an unreliable item. Therefore, RU cargo was considered 
part of the uncontrolled cargo flow and grouped with LRU 
data for the purpose of this work. Table VI, Appendix A, 
shows the number of RU shipments classified because of weight 
for each of the major PODs. 

It should be noted that no attempt was made to determine 
the number of shipments classified as RU for reasons other 
than weight (i.e., classified material). It was assumed 
that the volume of this cargo was of insignificant proportion 
when compared to the total yearly volume. 

In summary, the distribution audit provided a means of 
taking raw data and arranging it into a meaningful and use- 
ful format. It permitted the identification of the most 
active ports of debarkation which then became the focal 
point for further analysis. The grouping of the data into 
weekly intervals by POD indicated the tremendous volume 
fluctuations that occurred throughout the period and the ap- 
parent complexities that arise when attempting to forecast 
future volume flow. 

Up to this point, only the simple computer manipulation 
of data had been accomplished through the distribution audit. 
With this phase completed, there were still many questions 
to be answered. What type of data had been generated? Was 
there any correlation in the volume flow from week to week? 

The above questions are answered in the following section 
which begins the statistical analysis of the data. The ans- 
wers then provide statistical justification for applying 
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certain forecasting techniques which are covered in section 
C of this chapter. 

B. SERIAL CORRELATION TEST 

The flow of uncontrolled cargo volume into the Container 
Freight Division at MOTBA is time-series data. It is not 
uncommon for future observations in real time-series data to 
be dependent upon current observations. This dependence may 
be measured by analysis of the consecutive disturbances around 
a simple least-square-regression line. If these disturbances 
are found to be correlated from one observation to the next, 
and thus not random, serial correlation is said to be pre- 
sent. If the presence of serial correlation can be properly 
identified it can result in a more efficient estimation 
process in comparison to simple least -squares estimators. 

Such a process would give more weight to recent data for 
creating predictions of the values of future observations. 

Johnston [Ref. 9] lists three main consequences that 
result from applying straightforward least-squares formulas 
directly to observations that contain serially correlated 
disturbances. They are: 

1) The sampling variances of the intercept and the slope 
coefficient may be unnecessarily large in comparison to those 
obtained by other methods. 

2) The usual least-squares formulas for the sampling 
variances of the regression coefficients are no longer valid 
since a serious underestimate of the variance is likely to 
occur. 
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3) The sampling variances of the prediction will be 
unnecessarily large. 

1 . The Durbin-Watson d-Statistic 

When making forecasts using time-series data it 
would be a serious error to assume serial independence of 
the disturbance term. In order to avoid such an error, the 
Durbin-Watson d-Statistic [Ref. 6] provided a suitable test 
for the presence of serially correlated disturbances. The 
Durbin-Watson d-Statistic is defined by 

. tj2 (z t - z t-i j2 

where Z t (t=l , 2 , . . . ,n) denotes the residuals from a fitted- 
least-squares regression. To determine positive serial 
correlation the observed value of d is compared against 
lower and upper bounds of d^ and d^ that are tabulated for 
various values of n (the number of observations) and k (the 
number of explanatory variables) . The d-statistic for nega- 
tive serial correlation is calculated by taking 4 - d. 

This figure is also compared with d^ and dy Since the sign 
of the serial correlation was unknown, two-sided tests were 
conducted by combining single-tail tests. The test has four 
possible outcomes: (1) positive serial correlation of the 

data; (2) negative serial correlation of the data; (3) in- 
conclusive evidence of serial correlation with further 
observations being ideally required; or (4) no presence of 
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serial correlation (which implies that the disturbances are 
random) . 

The calculations were made using three different 
intervals. The first interval used weekly data resulting 
in 52 observations. The second interval used biweekly 
totals beginning with the observation for the first two 
weeks of FY '74. The total number of observations in this 
set was 26. Biweekly totals beginning with the observation 
of the second and third weeks of FY *74 (i.e., ignoring the 
first week's volume total) was used for the third interval. 
Thus there were only 25 observations resulting from this 
grouping scheme. 

The following example may help to clarify the above 
procedure. Given that the weekly volume input for POD, BA3 
for the first seven weeks of FY '74 is 10, 20, 30, 40, 50, 
60 and 70 cubic feet respectively, the first seven observa- 
tions for the weekly interval data set correspond according 
ly. The first three observations for the first biweekly 
data set will be 30, 70 and 110 cubic feet. However, the 
first three observations for the second biweekly data set 
will be 50, 90 and 130 cubic feet. 

The d-statistic test was implemented through Com- 
puter Program One. In addition to displaying the results 
of the calculations, the program displays the meaning of 
the d-statistic (i.e., positive, negative, inconclusive or 
none) . Computer Output One at the end of this thesis is a 
sample of the printout for Computer Program One. 
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2. Results of Serial Correlation Test 



While the Durbin-Watson d-Statistic was calculated 
for all 157 PODs , the d-s~tatistic for the 16 most active 
PODs was the focal point for the analysis. 

The results of the test for the weekly interval data 
set having 52 observations for each POD indicated positive 
serial correlation in the weekly volume inputs for nine PODs. 
The test was inclusive for four PODs and indicated no serial 
correlation for three others. There were no PODs with nega- 
tive serial correlation. In fact, there was no indication 
of negative serial correlation under any of the circumstances 
tested for all 157 PODs. The following table shows the test 



results using 


only the 


16 major 


PODs . The 


test was 


conducted 


at a level of 


significance of 0. 


05. 






Week per 


Week 




Results 




No 


Interval 


Started 


Positive 


Negative 


Incon 


Correlation 


1 


1 


9 


0 


4 


3 


2 


1 


5 


0 


5 


6 


2 


2 


7 


0 


5 


4 



Table VII gives the d-statistic results for each POD 
for the three intervals tested. 

Given the results of the d-statistic tests it was 
then decided to test for the aggregate presence of serial 
correlation in the data for the 16 major PODs. The null 
hypothesis for the test was 

V P 1 = p 2 " = p 16 = 0 
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where p^,...,p^ symbolically represents the serial correla- 
tion coefficients for the 16 major PODs . The alternate 
hypothesis was 

H^: not all p's = 0. 

This hypothesis is, in essence, testing the plausi- 
bility that the individual cases indicating positive serial 
correlation were Type I errors, i.e., a rejection of the 
null hypothesis, p^ = 0 , when it was true. This was formu- 
lated as a Bernoulli process with the probability of success 
(the indication of positive serial correlation) equal to the 
level of significance for the d-statistic test (0.05). 

Since a Bernoulli process was involved, the binomial 
distribution function was applied. Letting the probability 
of success, p, equal 0.05; the number of trials, n, equal 
to 16; and the specific number of successes, r, equal to 5, 

7 and 9 respectively, the decision rule was 

Accept H q if P(R>r|n,p) > 0.05 

Reject H q if P(R>r|n,p) < 0.05 

where R equals the unknown number of successes. 

Using cumulative probability distribution tables the 
following results were obtained: 

P (R>5 | 16,0.05) = 0.0009 
P (R> 7 | 16,0.05) = 0.0001 
P(R>9| 16,0.05) < 0.0001. 

Based on these results the null hypothesis was re- 
jected and it was assumed that significant serial correlation 
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surely existed somewhere in the data. Therefore, there was 
justification for using methods other than ordinary-least- 
squares regression to forecast the volume of cargo arriving 
at the Container Freight Division at MOTBA. 

C. FORECASTING METHODS 

The results of the Durbin-Watson d-Statistic indicated 
the presence of serial correlation and thus the serial de- 
pendence of the observations. It showed that the distur- 
bances around an ordinary-least-squares regression line 
would not fall equally on each side of the line. As the 
disturbances tended to fall more on one side of the line 
than the other, this implied that the most recent distur- 
bance was the best predictor of the subsequent distrubance. 
Also, by implication, the best predictor of the next obser- 
vation is the most recent observation or observations. This 
supplied justification for the application of forecasting 
methods other than ordinary- least- squares to place more 
emphasis on recent data rather than the minimization of the 
squared residual term. While there were a number of methods 
from which to choose, the methods used in this work were 
chosen because of their relevency to the problem and for 
their ease of calculation. 

This section looks at: the moving average method that 

is presently being used at MOTBA; the exponentially weighted 
average method; the exponentially weighted average method 
with adjustment for trend; and an adaptive response rate 
technique applied to the exponentially weighted average 
method . 
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These methods were evaluated using mean absolute devia- 
tion (MAD) and mean deviation squared (MDS) as the criteria 
for goodness. By computing the mean of the deviations it 
was possible to compare the various forecasting methods when 
they resulted in differing numbers of forecasts generated 
during the period. 

In addition, MAD and MDS were computed for one-, two- 
and three -week time horizons. For the one-week time horizon 
the forecast was subtracted from the actual input volume on 
a week by week basis. For the two- (three- )week time horizon 
the latest forecast was multiplied by two (three) and sub- 
tracted from the actual input volume for the two (three) 
weeks covered by the forecast. This was done on a week by 
week basis. 

The following example may help to clarify time horizon 
calculations. Say that the forecast for week 20 is 1,000 
cubic feet. If a three-week time horizon is being used, the 
forecast will be 3,000 cubic feet for weeks 20, 21 and 22. 

If the actual input volumes turn out to be 900; 1,100 and 
1,200 cubic feet respectively, the absolute deviation for 
the forecast is then 200 cubic feet (or 900 + 1,000 + 1,200 
- 3,000). Week 21 will result in a new forecast for weeks 
21, 22 and 23. The forecast for week 21 is multiplied by 
three and subtracted from the actual input volumes for weeks 
21, 22 and 23 to determine the absolute deviation for that 
particular forecast. Thus, for each forecast there is an 
absolute deviation that is computed. 
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What, then, was the purpose for using time horizons? 

One reason was the requirement for booking container space 
three weeks in advance and thus projecting a one-week fore- 
cast over a period greater than one week. One-, two- and 
three-week horizons were used to make comparisons of any 
effects that may be caused by changes in the booking re- 
quirements. It was also anticipated that an increase in 
the time horizon may have some effect on the results of the 
various forecasting methods due to a dampening effect on 
data with large fluctuations. This dampening effect is ex- 
plained later. 

It should also be noted that comparisons among the var- 
ious forecasting methods were only made within a given time 
horizon . 

Tables VIII and IX, Appendix A, show the MADs and MDSs 
generated by all of the forecasting methods used in this 
analysis . 

1 . Moving Average Method 

As previously explained, the CSS at MOTBA computes 
the previous three weeks average receipts in order to fore- 
cast the input volume that will be on hand in 21 days. The 
offerings for bookings that are made depend on three factors 
1) the volume of cargo on hand; 2) the volume of cargo 
booked; and 3) the forecasted volume. 

The three-week moving average is nothing more than 
a linear average. This average, multiplied by three is then 
used to forecast the volume receipts for the next three 
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weeks. In this method equal weight is given to each of the 
past three weekly volume totals. 

In order to evaluate the moving average method the 
MADs and MDSs were calculated for one-week to ten-week 
moving averages for the 16 PODs . In addition, the aggre- 
gate MAD and MDS for each POD were calculated for one- to 
ten-week moving averages. In other words, a specific moving 
average (from one to ten weeks) was applied to the volume 
inputs for each POD. The MAD and MDS for each POD were 
calculated and then summed to get the aggregate. The same 
calculations were made for two- and three-week time horizons. 

Computer Program Two was utilized to obtain the 
calculations for the analysis of the moving average method. 
Computer Output Two shows a portion of the output that was 
generated by the program. 

a. Data Analysis of the Moving Average Method 

In 13 of 16 cases a moving average other than 
the three-week moving average provided a better forecast 
when MAD was the criterion of goodness. When using MDS as 
the criterion, the result was 15 of 16 cases. The time in- 
tervals that provided the best forecast will be called the 
"optimal moving average." 

These evaluations were made using a one-week 
time horizon. Similar results were obtained using two- and 
three-week time horizons except that the optimal time inter- 
val for the moving average was not the same in each case. 

It was noted that in some cases the selection of 
the goodness criterion determined the optimal interval for 
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the moving average. This was the case for six PODs using 
a one-week time horizon and was true to a greater or lesser 
extent for the two- and three-week time horizons. 

The decrease in the MAD when a moving average 
other than a three-week moving average was optimal, varied 
from 4.5% to 18.5%. Similar reductions in MAD and MDS oc- 
curred when the other time horizons were used. 

The optimal time interval for a moving average 
that applied to all PODs was seven, ten and ten weeks for 
one-, two- and three-week time horizons, respectively, with 
MAD as the criterion. The same time intervals held true 
with MDS as the criterion. 

It was noted that the number of forecasts de- 
creased as the time interval for the moving averages in- 
creased. For example, in a three-week moving average, the 
first three weeks are averaged to make a forecast for week 
four. Thus, there were only 49 forecasts made for the 
period. For a ten-week moving average there were only 42 
forecasts. The use of two- and three-week time horizons 
resulted in even fewer forecasts. Comparison among the 
various moving averages for a specific time horizon was 
possible since mean absolute deviation and mean deviation 
squared were used as the criteria for goodness. 

As the time interval for the moving average in- 
creased, the greater was the chance for large fluctuations 
in the data to be dampened by the initial forecast. For 
example, the volume inputs for the first six weeks of FY 
*74 for TA1 were 2,611; 4,307; 1,079; 3,360; 6,765 and 
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2,429 cubic feet, respectively. These fluctuations were 
dampened when moving averages greater than three weeks were 
applied to the data (with a five-week average being opti- 
mal) . This dampening effect stressed the fact that a three- 
week moving average placed more emphasis on recent data. 
Longer time intervals placed less emphasis on recent data 
and considered older data as well. 

It was also noted that the MAD and MDS decreased 
as the time horizon increased. The reason was that the 
two- and three-week time horizons also tended to cancel 
the volume fluctuations. A look, again, at Figure 1 shows 
the very large volume fluctuations that occur in the POD 
data and how it is possible for the dampening effect to take 
place . 

2 . Exponentially Weighted Average Method 

In lieu of the results of the serial correlation 
test, a viable alternative to the moving average method is 
the exponentially weighted average (EWA) . In its simplest 
form an EWA is based on a period by period adjustment of the 
most recent forecasted average by adding or subtracting a 
fraction of the difference between the actual volume in the 
current period and the latest forecasted average volume. 

It differs from the moving average in that the weight given 
the most recent observation can be varied while also giving 
consideration to all past data. EWA is represented by 

F t - a V t ♦ (1 - a) Vi 
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where is the forecast average, V t is the actual input 
volume and a is the exponential smoothing constant which 

must be between zero and one. The forecast for the next 

* 

period F t+ ^ can be taken directly from the computed value 
of F^. This can be justified because trends and seasonal 
adjustments are not accounted for in the model [Ref. 7], 

The exponentially weighted average method is also 
known as the exponential smoothing method. Both names are 
used in this thesis. 

The calculations for the EWA method were implemented 
by Computer Program Three. The MAD and MDS for each POD 
were calculated for a-levels between -0.10 and +1.00 in in- 
crements of 0.01. I'n addition the MAD and MDS were computed 
for a-levels that were applied to all 16 PODs . In other 
words a specific a-level was applied to data for each POD. 
The MAD and MDS for each POD were calculated and then summed 
to obtain an aggregate total. The a-levels which generated 
the smallest MAD and MDS are referred to as the "optimal 
a-levels for aggregate data." It should be noted that the 
optimal a-levels for aggregate data may differ depending on 
the criterion of goodness. 

A canned computer subroutine, SUBROUTINE OSPLOT, 
available at the computer center, was used to plot the 
a-levels versus the MAD for each POD. The plot displayed 
how the MAD reacted to various levels of a. Thus it was 
possible to visually determine the optimal a-level for 
aggregate data with MAD as the criterion of goodness. 
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Computer Output Three -Two shows the MADs and MDSs 
calculated for various a-levels. Computer Output Three-One 
is a sample of the OSPLOT for four of the PODs . 
a. Data Analysis of EWA Method 

With MAD as the criterion for goodness, the EWA 
method resulted in a decrease in MAD in 14 of 16 cases when 
compared to the three -week moving average. The decrease 
ranged from 4.9% for RJ3 to 12.4% for XJ1. However, when 
EWA was compared to the optimal moving averages for the 
individual PODs it produced a smaller MAD in only 4 of 16 
cases. The a-levels varied from 0.04 to 0.60 for the 16 
PODs. 

The optimal a-level for aggregate MAD data was 
0.36. The aggregate MAD for the EWA method was 59,387 
cubic feet as compared to 60,172 cubic feet for the three- 
week moving average and 59,209 cubic feet for the optimal 
moving average. 

The preceding analysis was made from calcula- 
tions using a one-week time horizon. Similar results were 
obtained using two- and three-week horizons. It is interest- 
ing to note that there were changes in the a-levels for most 
PODs depending on the time horizon. One POD (UM4) produced 
negative a-levels (-0.02 and -0.01) for two- and three- 
week time horizons, respectively. The optimal a-levels for 
aggregate data were 0.35 for the two-week horizon and 0.33 
for the three-week horizon. As was the case with the moving 
average method, the aggregate MADs decreased as the time 
horizons increased due to the dampening of the fluctuations. 
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With MDS as the criterion for goodness the EWA 
method resulted in a decrease in MDS in 12 of 16 cases when 
compared to the three -week moving averages. When compared 
to the optimal moving average, there was a decrease in only 
5 of 16 cases. In each case the optimal a-level for MDS 
was different from the optimal level for MAD; although the 
optimal a-level for aggregate data was approximately equal 
for the MDS (0.37) and the MAD (0.36). The aggregate MDS 
for the EWA method was less than the MDS for the three-week 
moving average but greater than the MDS for the optimal 
moving average. 

A one-week time horizon was used in the above 
analysis for MDS. As was the case with the MAD analysis 
of individual PODs , two- and three-week time horizons 
generated similar results. However, for a two-week time 
horizon of aggregate data, EWA generated a slightly smaller 
MDS than that generated by the optimal moving average. 

This was not the case for the three-week time horizon, 
however . 

The way in which the EWA method was applied re- 
sulted in no forecast deviations being computed until the 
third week. This resulted in some smoothing of the data but 
not as much smoothing as occurred with moving averages with 
intervals greater than three weeks. When applying the EWA 
method to data with large fluctuations at the beginning of 
the year, there was a lag in the forecast until these fluc- 
tuations were dampened over time. This was the apparent 
reason for the seven-week moving average being a better 
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forecaster than the EWA for a one-week time horizon and a 
ten-week moving average being better for a two- and three- 
week time horizon. 

In summary, the EWA method provided a means of 
emphasizing past data as well as placing some significance 
on the most recent data by the use of the smoothing con- 
stant. The EWA method provided a better forecast than the 
three-week moving average; however, the decrease in the MAD 
for aggregate data was only 2.21. For individual PODs this 
decrease was as high as 18 %, but for other PODs there was 
an increase in the MAD when EWA was used. 

3. Exponentially Weighted Average Method with Trend 

Adj ustment 

Neither the moving average or the EWA took into 
consideration any adjustment for trends that may exist in 
the data. Also, the slopes of the least-squares -regression 
line for some PODs used in the d-statistic calculations gave 
an indication of the possibility of trend effects in the 
data. Therefore, the EWA was modified to take into account 
the adjustment for trend. 

The difference between the volume forecast averages 
from week to week is the apparent trend. Exponential smooth- 
ing can be applied to this trend just as it was applied to 
the forecast average. This trend adjustment is represented 
by 

T t - “< F t - F t-P + O - a) T t . r 

By combining this trend adjustment with the new forecast 
average, the forecast for the next period is then 
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or 



F* 

t+1 



= V. 



(1 - a) F 



t-1 



— T 
a t-1 




where F t+ ^ is the forecast for the exponentially weighted 
average with trend adjustment (EWA-TA) [Ref. 7], 

Computer Program Four was used to make the statisti- 
cal calculations for the EWA-TA. Except for the trend ad- 
justment, the computational procedures were the same as those 
for the EWA method. Once again, SUBROUTINE OSPLOT was used 
to plot a-levels versus MAD for each POD. 
a. Data Analysis of EWA-TA 

Using a one-week time horizon and MAD as the 
criterion, a decrease in MAD resulted in only 3 of 16 cases 
when comparing EWA-TA to EWA. Only in the case of RG1 was 
there any significant change (an 8% decrease). However, 
this POD was an exception to the rule of fluctuating volume 
inputs. Its weekly volume inputs toward the last half of 
the year had remained within a relatively small range. In 
only one case did EWA-TA provide a better forecast than any 
of the methods previously discussed. However, the improve- 
ment was minimal. In half the cases the difference in MADs 
between EWA-TA and EWA was less than 50 cubic feet. Also, 
the a-level dropped significantly in all 16 cases. The op- 
timal a-level for aggregate data was 0.16; however, the MAD 
was greater by 1,904 cubic feet. 
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With MDS as the criterion for goodness the re- 
sults were similar. Improvement on the EWA method occurred 
in only two cases. In most cases the a-levels were differ- 
ent from those with MAD; however, the optimal a-level for 
aggregate data changed by only 0.01 to 0.17. The aggregate 
MDS was the largest of all the methods applied thus far. 

This was true for all time horizons. 

Very few changes occurred when two- and three- 
week time horizons were applied. Once again, UM4 produced 
negative a-levels. The optimum a-levels for aggregate data 
did not change significantly. The aggregate MAD and MDS 
decreased in the same proportion as the aggregate totals for 
the other methods. 

In summary, the exponentially weighted average 
method with trend adjustments did not provide an improved 
forecast over the moving average method or the exponentially 
weighted method without trend adjustments. Apparently there 
was no definite trend in the volume inputs for the PODs . It 
is doubtful that the method can be effective when applied 
to data covering only a one year period. 

4 . Exponentially Weighted Average with an Adaptive 
Response Rate 

Because of the large volume fluctuations in the data, 
the application of an adaptive response rate (ARR) provided 
a plausible alternative to a simple exponentially weighted 
average. The ARR provides an automatic means of detecting 
sudden volume changes through the use of a tracking signal. 
Trigg and Leach [Ref. 18] define the tracking signal as 
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Tracking Signal 



Smoothed Error 

Smoothed Absolute Error 



where the error is the difference between the forecasted 
volume and the actual volume inputs in each period. The 
error was smoothed by choosing a value for gamma, y, between 
-1 and +1 such that 



New Smoothed Error = (1-y) old smoothed error + y error, 



New Smoothed Absolute 
Error 



(1-y) old smoothed absolute error 
+ y absolute error. 



The tracking signal was then applied to the EWA methods with- 
out trend adjustment by letting the smoothing constant equal 
the absolute value of the tracking signal. Thus, it ap- 
peared that ARR would provide a rapid response to large 
volume changes whereas the constant model would lag behind 
these sudden changes. 

The adaptive response rate was also applied to the 
exponentially weighted average with trend adjustment. It 
resulted in very large deviations because of the conflict 
between the trend adjustment and the response rate. They 
were apparently working against each other; and therefore, 
this method will be ignored in this paper. 

MADs and MDSs for the 16 major PODs were calculated 
for y-levels between -1.00 and +1.00 in increments of 0.01. 
Optimal y-levels for aggregate data were also computed. 

These calculations were implemented by Computer Program 
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Five. Plots of ylevels versus MAD were made through the 
use of the special subroutine 

a. Data Analysis of EWA with ARR 

With a one-week time horizon and MAD as the 
criterion, the ARR resulted in the smallest MAD in 6 of 16 
cases when compared to all other methods previously dis- 
cussed. The reduction in MAD was less than 51 in five of 
these cases. The optimal y-level for aggregate data was 
-0.05. The aggregate MAD at this level was smaller than 
the MADs for the three-week moving average and the EWA-TA; 
but it was larger than the MADs for the optimal moving 
average and the EWA method. 

A two-week time horizon produced similar re- 
sults on the individual POD as well as aggregate level. 
Individual POD results were similar with a three-week time 
horizon. For aggregate data only the optimal moving aver- 
age resulted in a smaller MAD. However, the reduction was 
less than 1%, and for all practical purposes it was equal. 

With a one -week time horizon and MDS as the 
criterion, the ARR produced a lower MDS in 5 of 16 PODs 
when compared to all other methods. Two- and three-week 
time horizons had similar results on the POD level. For 
aggregate data, the ARR was unable to generate an MDS smaller 
than any of the other methods for the three time horizons. 

In summary, the ARR provided a more complicated 
approach to exponential smoothing than any of the other 
methods. On the individual POD level it provided a better 
forecast in some cases; yet, the improvement was minimal. 
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The time horizon and the criterion for goodness were signif- 
icant factors in determining whether the ARR generated the 
better or best forecast. 

While the purpose of the ARR was to respond more 
rapidly to large fluctuations in the data, it would be more 
appropriate to state that it was designed to respond to 
large fluctuations which resulted in a shift in the data. 
While there were large fluctuations in the data there were 
no apparent shifts in the data. Plots of the weekly volume 
inputs indicated that there were large fluctuations in the 
data, but that these fluctuations did not result in any 
significant shifts in the data. 

5 . Summary of Forecasting Methods 

Four forecasting techniques were applied to the 
data: 1) the moving average method; 2) the exponentially 

weighted average method; 3) the exponentially weighted 
average method with trend adjustment; and 4) the exponen- 
tially weighted average method with an adaptive response 
rate. Past data was utilized to forecast future inputs in 
each of the methods. Each method provided a means of 
varying the amount of emphasis placed on the most recent 
data, forecast average, or error. 

In order to evaluate the forecasting methods, the 
mean absolute deviation and mean deviation squared were 
calculated and used as the criterion of goodness. The mean 
of the deviations was used since each method resulted in 
a different number of forecasts being made for each period. 
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The MADs and MDSs that were generated by the various 
forecasting methods indicated the uniqueness of the volume 
inputs for each POD. The method providing the best fore- 
cast for a POD depended upon the criterion of goodness and 
the time horizon. 

Tables VIII, IX, X and XI, Appendix A, are a pre- 
sentation of the data that has been previously discussed. 
Table VIII gives the MAD of the various methods for the 16 
major PODs in one-, two- and three -week time horizons. 

Table IX gives the MDS in the same format. Tables X and XI 
present the MAD and MDS for aggregate totals respectively. 

From the data in these tables it can be seen that 
lower MADs and MDSs may result from increasing the time 
horizon. It has been suggested that the lower deviations 
were a mathematical result of the two- and three-week time 
horizons averaging the volume fluctuations. The reduction 
continued for each increase in the time horizon in 8 of 16 
PODs. In only one case did a decrease not result from an 
increase in the time horizon from one to two or three weeks. 
In other words, in 15 of 16 cases there was a reduction in 
MAD when using a two- or three-week horizon. 

The moving average method provided a means of damp- 
ening the initial large volume fluctuations in the data. 

Large volume fluctuations during the first quarter of FY 
'74 could be dampened by using a 10-week moving average. 

Thus, the first deviation was not computed until the eleventh 
week . 
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The other methods that were applied to the data be- 
gan computing deviations as early as the third week. With 
large volume fluctuations at the beginning of the period it 
was possible that there was a lag in the forecast which 
initially resulted in large deviations. The plots of the 
weekly inputs showed that 8 of 16 PODs had large volume 
fluctuations at the beginning of the period. • 

An overall view of the forecasting methods indicated 
that no one particular method dominated the others in terms 
of providing the best forecast on a POD or an aggregate 
level. The difference between the largest and smallest MAD 
was less than 100 cubic feet for one POD (UL7) . For another 
POD (UM4) three different methods were optimal depending on 
the time horizon. These and other factors demonstrated the 
random flow of cargo arriving at the CSS. 

6 . Optimal Forecasting Method Selection 

The results of the various forecasting techniques 
applied to the data provide a basis for the selection of an 
optimal method. Through the application of these techniques 
the following alternatives were available for forecasting 
the weekly volume inputs of the PODs : ^ 

(1) a three-week moving average, 

(2) a moving average with an optimal time interval 
applied to all PODs. 



■*■566 Tables X and XI, 



Appendix A. 
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(3) a moving average with the optimal interval for 
each POD applied to its weekly inputs (i.e., the time in- 
tervals vary among the PODs) , 

(4) an exponentially weighted average with the op- 
timal a-level applied to all PODs, 

(5) an exponentially weighted average with the op- 
timal a-level for each POD applied to its weekly inputs 
(i.e., the a-level varies among the PODs), 

(6) an exponentially weighted average with trend 
adjustment with the optimal a-level applied to all PODs, 

(7) an exponentially weighted average with trend 
adjustments with the optimal a-level for each POD applied 
to its weekly inputs, 

(8) an adaptive response rate with the optimal 
y-level applied to all PODs, 

(9) an adaptive response rate with the optimal 
y-level for each POD applied to its weekly inputs. 

The selection of the optimal method required an 
analysis of the methods used to calculate the goodness 
criterion and the application of Student's t-test in order 
to determine if the -difference of the deviations among the 
methods were significantly different. 

The results of the one to ten-week moving average 
indicated that a time interval greater than three weeks 
generated smaller deviations than the three-week moving 
average for 13 of 16 PODs. This ratio held true for both 
criterions and for all time horizons. Thus, if a moving 
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average is to be applied to the data, a time interval greater 
than three weeks is preferred. 

There is one shortcoming in comparing the moving 
average method with the others. The first deviation for a 
moving average was not computed until the week following 
the time interval of the average while the deviations were 
computed beginning the third week for the other methods. 

Thus, the longer time intervals tended to have a greater 
initial smoothing of the forecasts while there was a lag 
in the smoothing for the other methods. This fact should 
be kept in mind when comparing the aggregate totals in 
Table X, Appendix A. Because of the manner in which the 
deviations were computed no comparisons were made between 
the optimal moving average and the other methods. 

The exponentially weighted average offers an alter- 
native to a moving average. Firstly, it did generate smaller 
MADs and MDSs than a three-week moving average. A simple 
t-test using MAD as the criterion and a one-week time hori- 
zon indicated that alternatives (4) and (5) are significantly 
different from alternative (1) at the 0.05 confidence level; 
and, therefore, (4) and (5) are preferred to (1). At the 
same confidence level, the t-test indicated that (4) and 
(5) are different; and, therefore, alternative (5) is pre- 
ferred to (4) . Not until a confidence level of 0.005 is 
used, can it be said that alternatives (4) and (5) are not 
significantly different. 

The exponentially weighted average method with trend 
adjustment (alternatives (6) and (7)) generated aggregate 
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MADs and MDSs that were larger than those generated by the 
EWA without the trend adjustment. Based on this fact, al- 
ternatives (4) and (5) are preferred to alternatives (6) 
and (7) . 

The exponentially weighted average with an adaptive 
response rate (alternatives (8) and (9)) provided another 
alternative to the moving average method. For 6 of 16 PODs 
the ARR generated the smallest MAD for one- and two-week 
time horizons. On the aggregate level alternatives (8) and 
(9) generated slightly smaller MADs than did alternatives 
(4) and (5). A t-test indicated that alternatives (8) and 
(9) are different; and, therefore, (9) is preferred to (8). 
The test was conducted at the 0.05 confidence level. 

With alternatives (5) and (9) generating smaller 
MADs and MDSs than the three-week moving average, the ques- 
tion was which alternative is preferred? Again, a simple 
t-test at the 0.05 level of confidence was applied. At this 
level, the test indicated that alternatives (5) and (9) are 
not significantly different. With this fact in mind, al- 
ternative (5) would be preferred to (9) because of compu- 
2 

tational ease. 

While it was possible to use a simple t-test to ar- 
rive at a preferred method of forecasting, the results of 
the various methods have indicated the true randomness of 
the cargo flow into the Container Stuffing Station. No one 



2 

Table XII, Appendix A, shows a sample of the t-test 
calculations . 
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particular method dominates the others in its ability to 
accurately forecast the weekly volume inputs for all of the 
major PODs . The satisfaction of arriving at a preferred 
method must come from the fact that this method does result 
in deviations that are smaller than those generated by a 
three -week moving average and is computationally easier and 
less sensitive to the data than another method that is com- 
parable in its forecasting accuracy. 

Additionally, the minimum MADs achieved using alter 
native (9) occurred at y-levels ranging from -1.00 to +0.49 
Also, the data for most PODs possessed several maxima and 
minima of MAD as y-levels were varied. This revealed the 
highly sensitive nature of the results of forecasts to the 
input data when alternative (9) was used. Should the char 
acteristics of these data change in future observations, 
errors resulting from the continued use of previously opti- 
mal y-levels could become substantial. This was contrasted 
to the more narrow range of a-levels (from 0.02 to +0.77) 
for alternative (5) and the achievement of a single minimum 
MAD for most cases. 

Computer Output Three -One and Computer Output Five- 
One are representative of these relationships. 

7 . Forecast Variance Estimation 

To this point, this thesis has been concerned with 
the choice of an optimal method of forecasting the volume 
inputs to the CSS, using mean absolute deviation and mean 
deviation squared as the criterion of goodness. The appli- 
cation of any method chosen requires that the statistical 
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variance of the forecast be known or estimated such that 
confidence in the forecast may be assessed. The knowledge 
of the variance can assist the booking of vans by providing 
a means of assessing the impacts and costs of conservative 
versus aggressive booking policies and the probable effects 
on average age of cargo, cancellations, cargo volume on the 
floor, cube utilization and other policies. Just as the 
expected volume was forecast by some method, so must the 
expected variance be forecasted or estimated. Two possible 
procedures for arriving at variances for the forecasts are 
suggested. 

The first possible method of estimating the variance 
and standard deviation of volume inputs for a POD could be 
to project forward historical information for one period 
to the next. If this method were chosen, the length of the 
period would dictate how much averaging out of fluctuations 
would be accomplished. Conversely, should the fluctuations 
be tied to some major policy decisions occurring during a 
period, variance estimates could be seriously in error if 
the period chosen was too long. It is suggested that per- 
haps a period of six months to one year, providing 26 to 
52 weekly deviation observations. Thus, the weekly volume 
input variance for a POD calculated for one period could be 
applied to all weekly forecasts for that POD in the ensuing 
period . 

Utilization of a smoothed square deviation computed 
simultaneously with each volume forecast is another possible 
method of estimating the variance. As the variance of each 
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forecast is related to the square of the deviation between 
the observed and forecasted volume for the period, this 
value could be smoothed from one forecasting period to the 
next using the simple exponentially weighted average tech- 
nique. The estimated standard deviation for the next fore- 
cast could be calculated as the square root of the smoothed 
squared deviation. The value of the smoothing constant 
(alpha) for this calculation could be selected depending on 
the weight to be placed on the most recent data (usually 
between 0.1 and 0.3 according to Brown [Ref. 3]). Intui- 
tively, the selection of a larger alpha could be tied to 
more aggressive booking policies by trying to detect recent 
trends while smaller values would provide more conservative 
long-run averages. In any case, this technique provides 
continual updating of variance estimates. 

Neither of these methods was tested persuant to 
this thesis due to the limitation of having only one year's 
data and the lack of any apparent measure of goodness that 
would indicate the relative optimality of one method over 
the other. 

Whether or not either of these methods is used, 
some estimates of the variance associated with a forecast 
prediction is imperative to establish consistency between 
the booking policy and other policies concerned with the 
movement through and storage within the Container Stuffing 
Station . 



74 



D. CRITIQUE OF PROCEDURES AND DATA ANALYSIS 

The analysis leading up to and including the forecasting 
techniques have been presented in detail in the preceding 
sections of this chapter. This section will evaluate the 
procedures used to make the analysis. 

1 . Distribution Audit 

The distribution audit was a necessary first step 
in the data analysis. Although there was a multitude of 
information available concerning the cargo arriving at MOTBA, 
only the number of shipments, their volume and their weight 
was of concern to this work. The audit took the raw data 
and placed it in a meaningful format. 

The distribution audit helped to establish the fact 
that the Vietnam conflict generated a tremendous cargo flow. 
Therefore, it was decided not to use the six months of data 
from FY '73. This left 52 weeks of data for FY *74, pro- 
viding a reasonable number of observations. Another alter- 
native would have been to exclude the data for the first 
two months of FY '74, because this data had large volume 
fluctuations on a scale similar to FY '73 data. The ex- 
clusion of this data may have resulted in the generation of 
lower deviations for any or all of the forecasting techni- 
ques that were applied. 

Since there was a requirement to book container 
space aboard vessels three weeks in advance, the data was 
grouped in weekly intervals resulting in 52 observations 
for the period. Although biweekly data could have been used, 
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it was more convenient to use weekly intervals for the data 
analysis. Intervals greater than two weeks would have re- 
sulted in too few observations. 

The distribution audit also identified the number 
of shipments classified as release unit because of weight. 
These shipments were considered part of the random volume 
flow based on the fact that they comprised a very small 
percentage of the total number of shipments. 

In short, the distribution audit provided a start- 
ing point for the procedures that followed. It simplified 
the data by placing it in meaningful perspective. 

2 . Serial Correlation Test 

The data was analyzed to reveal the presence of 
serial correlation. The Durbin-Vvatson d-Statistic Test is 
commonly used for the purpose. The calculations for the 
test were easily accomplished by computer programming. The 
data was tested using weekly and biweekly intervals. Tri- 
weekly intervals were not used since they would have re- 
sulted in only 17 observations. The results of the test 
justified the application of the forecasting methods pre- 
viously discussed. 

3 . Forecasting Methods 

a. Moving Average Method 

The moving average was the most obvious method 
to initially apply to the data. It provided a simple means 
of forecasting future volume inputs. Less weight was given 
to the weekly input as the time interval for the moving 
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average was increased. Time intervals from one to ten weeks 
were evaluated. To carry the evaluation beyond ten weeks 
would have been superfluous since time intervals beyond ten 
weeks result in minimal changes in the forecasts for each 
new input of data. As the moving average intervals become 
large, the process approaches the use of the mean of the 
volume inputs over a certain time period (for example, a 
quarter) as a forecast for the next similar period. The 
calculations indicated that increased smoothing of the data 
was accomplished by increasing the time interval of the 
moving average beyond three weeks. 

b. Exponentially Weighted Average Methods 

While previous discussions may have indicated 
that there were four forecasting techniques applied to the 
data, in reality there were only two -- the moving average 
and the EWA plus two variations. 

The EWA method, unlike the moving average, con- 
siders all past data to some extent. Through the applica- 
tion of the smoothing constant, the emphasis placed on the 
most recent data can be varied. 

One variation of the EWA was to adjust for trend 
effects in the data. The apparent trend was then smoothed 
and combined with a smoothed forecasted average to make a 
forecast for the following week. 

The second variation was the application of an 
adaptive response rate to the EWA. This particular ARR used 
a smoothed error term to calculate a tracking signal. Other 
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methods for calculating a tracking signal may be available; 
but they were not investigated. 

When the EWA methods were applied to the data, 
deviations were calculated beginning with the third week. 

For the moving average method, deviations were not computed 
until the week following the length of the time interval. 

In order tQ evaluate the methods it may have been more ap- 
propriate to begin calculating the deviations for the EWA 
methods in the eighth, ninth or tenth week. 

4 . Data 

The time period covered by the data was another 
limitation in the application of the forecasting techniques. 
In general, time-series analysis covers a three to four year 
period with data being grouped by months or even quarters. 
Conditions for this analysis were far from being optimal. 

The fact that there was only a year's data also pre- 
cluded investigation of seasonal or cyclic effects. Again, 
a three to four year period is usually required before 
smoothing techniques can be applied to the data. 

In summary, the breakdown of the yearly data into 
weekly intervals constituted the most logical approach to 
the problem. While a year of data did not represent a long 
period of time, it did represent the best that was available 
at the time. 
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V. CONCLUSION 



This thesis concerned itself with the forecast of the 
volume of cargo flowing into the Container Stuffing Station 
at MOTBA destined for a particular port of debarkation. 

The data used in this work was taken from actual operations 
at the CSS during FY '74. Although every system is capable 
of generating errors, it was assumed that any errors were 
minimal and that the data was correct. 

Prior to utilizing any forecasting techniques, the data 
was analyzed in three phases: a distribution audit; the 

grouping of data into weekly intervals by POD; and a test 
for serial correlation. 

This thesis made use of the distribution audit of 
McCarthy and Carter. The purpose of the audit was to sim- 
plify the large volume of raw data generated by the opera- 
tions at MOTBA. The magnitude of the operations at the 
CSS is illustrated by the following information. During 
FY '74 there were a total of 124,592 shipments weighing 
over 176 million pounds and with a total volume of 10,070,724 
cubic feet. These shipments were destined for 3,356 consig- 
nees in 157 overseas ports. Sixteen major ports accounted 
for 931 of the shipments, 93% of the volume and 91% of the 
weight. These 16 PODs then became the focal point for 
further analysis. 

The grouping of the data into weekly intervals by POD 
indicated the large volume fluctuations that occurred for 
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most of the 16 PODs under consideration. Plots of the 
weekly volume inputs revealed the lack of any set pattern 
of cargo flow from POD to POD. Before making any decisions 
as to which forecasting techniques to apply, it was neces- 
sary to determine whether the data was random in nature or 
correlated from one observation to the next; that is, 
serially correlated. The Durbin-Watson d-Statistic Test 
satisifed this requirement. A computer program was written 
to compute the d-statistic and to compare it with critical 
values. Using weekly volume inputs the test indicated that 
9 of 16 PODs had serial correlation in the data. Another 
test was performed by grouping the data into biweekly in- 
tervals. Similar results were obtained. The preceding 
results were then tested to determine if there was serial 
correlation present in the data. The test was affirmative. 

The results of the serial correlation test supported 
the use of forecasting methods other than ordinary least- 
squares trending. The simplest and most logical forecast- 
ing methods to apply to the data were the moving average 
method and the exponentially weighted average method. Two 
variations were applied to the exponentially weighted aver- 
age resulting in a total of four methods. 

In order to evaluate the forecasting methods on an 
equal basis, the mean absolute deviation and the mean de- 
viation squared were used, the deviation being the differ- 
ence between the forecast and the actual input. Four 
computer programs were written to calculate the MADs and 
MDSs for the various forecasting methods. 
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Moving averages from one to ten weeks were applied to 
the data. The MAD and MDS for each time interval were 
calculated by computer program. The EWA method used a 
smoothing constant to generate the forecast. Various levels 
of the constant were applied to the data. A variation of 
the EWA calculated the apparent trend of the data. This 
apparent trend and the forecasted average were then smoothed 
and added to get the forecast for the following week. An- 
other variation to the EWA method was the application of an 
adaptive response rate. The ARR uses a tracking signal 
which is computed by taking the sum of the differences be- 
tween the forecast and the actual input volume for each 
week and dividing this total into the sum of the absolute 
values of those differences. As each error becomes avail- 
able it is smoothed along with the old error. The smooth- 
ing constant to be used in forecasting is then set equal to 
the absolute value of the tracking signal and then the EWA 
is applied as before. 

The deviation statistics resulting from the application 
of the above forecasting methods to the actual data with a 
one-week time horizon were tested using the Student t-dis- 
tribution. This test revealed that the EWA and ARR methods 
using the individual POD's optimal smoothing constants were 
statistically the same and better than all other methods. 

The EWA technique was then chosen as preferred based on 
its computational ease relative to ARR. 

The selection of the EWA method offered some advantages 
over the three-week moving average being used at the CFD. 
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First, EWA was computationally no more difficult. It de- 
rived its advantage by responding to changes in volume more 
uniformly than the abrupt alterations that the three-week 
moving average caused in the forecasts. Also, the accumula- 
tion of a historical data base in the smoothed forecast 
average increased the accuracy of the EWA method as the time 
horizon increased. Decreases of 4.91, 8.0% and 10.0% in 
mean absolute deviations for the one-, two-, and three-week 
time horizons using the EWA method evidenced the increase 
in accuracy. For this application and data, there can be 
little doubt that the use of the EWA provided improvement 
in the forecasting techniques. 

The other forecasting methods retained their relative 
merits. Greater accuracy was achieved using moving averages 
by increasing the lengths of the averaging interval at the 
expense of increasing the computational difficulty. The 
absence of any apparent trend in the data resulted in lack 
of any additional benefits from trend adjustments. While 
the adaptive response rate technique showed no real advan- 
tage over the EWA method for this data, its ability to 
respond rapidly to major policy changes makes it a viable 
alternative at the expense of increased computational dif- 
ficulty. As more data becomes available, all of these 
methods deserve additional consideration and thoughtful 
attention as the results of this work are based on only a 
year's data. 

Although significant increases in forecasting have been 
achieved, there remains in the data large amounts of randomness 
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causing significant errors. Its random nature, without 
trends or significant shifts* have l.ed to the evolution of 
the exponentially weighted average method as the preferred 
choice. As more data is accumulated the other methods may 
be reexamined. The optimal smoothing constants should be 
updated periodically to account for changes in the data and 
promote long run efficiency of the method. However, the 
most important safeguard against costly mistakes due to 
random effects in the input volume is careful managerial 
scrutiny of each forecast generated and used to make a 
booking. Judgment on the reasonableness of a prediction 
makes the manager the final and ultimate smoothing factor 
in the utilization of the forecast. 
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APPENDIX A 



WEEK 




SHIPMENTS 




VOLUME 


WEIGHT 


27 




3,208 




220,991 


3,925,310 


28 




5,091 




413,602 


7,164,543 


29 




4,598 




336,189 


6,062,987 


30 




4,884 




368,350 


6,038,153 


31 




5,707 




344,109 


6,312,502 


32 




4,727 




372,533 


6,281,837 


33 




4,726 




358,880 


6,640,755 


34 




4,214 




303,929 


5,393,130 


35 




5,681 




378,598 


7,275,388 


36 




5,031 




324,176 


6,206,201 


37 




5,045 




348,497 


6,760,475 


38 




4,638 




301,952 


5,808,959 


39 




4,384 




303,536 


5,719,112 


40 




4,047 




300,089 


5,138,948 


41 




4,646 




335,702 


6,322,970 


42 




4,753 




343,370 


6,425,117 


43 




4,317 




314,323 


5,969,897 


44 




3,935 




281,744 


5,190,323 


45 




3,906 




253,313 


4,840,206 


46 




3,318 




248,060 


4,473,395 


47 




3,660 




288,476 


5,085.827 


48 




2,957 




231,423 


4,205,199 


49 




3,498 




320,646 


5,902,697 


50 




3,338 




285,413 


1,150,959 


51 




3,675 




387,533 


5,980,555 


52 




3,579 




433,497 


6,248,080 


Note : 


(1) 


Volume is 


in cubic 


feet 






(2) 


Weight is 


in pounds 




Table ] 


Gross Weekly 


Inputs . 


Weeks 27 


through 52, FY ' 
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WEEK 


SHIPMENTS 


VOLUME 


WEIGHT 


1 


2,444 


261,448 


3,921,219 


2 


2,797 


287,444 


4,527,403 


3 


3,141 


278,991 


4,997,954 


4 


2,758 


232,296 


3,571,542 


5 


3,414 


284,814 


5,151,369 


6 


4,044 


308,446 


5,844,244 


7 


965 


88,024 


1,484,204 


8 


2,905 


208,179 


3,908,164 


9 


2,641 


196,694 


3,625,615 


10 


1,994 


164,170 


3,017,826 


11 


2,425 


206,833 


3,837,876 


12 


2,787 


215,422 


4,308,528 


13 


2,485 


199,702 


4,533,820 


14 


2,541 


192,863 


3,600,097 


15 


1,784 


135,561 


2,378,525 


16 


2,653 


194,927 


3,742,275 


17 


1,983 


150,028 


2,857,312 


18 


1,387 


111,154 


2,199,669 


19 


2,340 


208,560 


3,761,167 


20 


2,411 


184,520 


3,200,648 


21 


1,361 


105,629 


1,965,960 


22 


2,414 


149,722 


3,663,732 


23 


2,927 


203,497 


3,864,209 


24 


2,754 


199,888 


3,695,825 


25 


2,757 


189,969 


3,657,957 


26 


1,523 


130,513 


2,592,805 



Note: (1) Volume is in cubic feet 

(2) Weight is in pounds 



Table II. Gross Weekly Inputs. Weeks 1 through 26, FY '74. 
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WEEK SHIPMENTS 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

Note : 



VOLUME 


WEIGHT 


135,843 


2,186,285 


232,446 


3,841,574 


260,123 


4,468,724 


179,400 


3,078,738 


241,953 


3,793,319 


242,661 


4,034,761 


211,740 


3,470,228 


139,665 


2,459,631 


267,702 


4,907,835 


164,209 


2,725,845 


143,360 


2,385,097 


145,621 


2,851,101 


172,236 


3,151,837 


170,630 


2,921, 166 


201,845 


3,388,849 


195,070 


3,211,454 


191,623 


3,135,205 


198,084 


3,052.520 


204,956 


3,269,329 


193,582 


3,264,435 


177,121 


3,355,770 


96,791 


1,628,575 


185,867 


2,771,581 


194,153 


2,716,456 


229,928 


3,602,602 


188,520 


2,631, 603 



(1) Volume is in cubic feet 

(2) Weight is in pounds 



1,530 

2,956 

3,058 

2,338 

2,580 

2,814 

2,667 

1,585 

3,007 

2,083 

1,994 

2,198 

2,557 

2,247 

2,060 

2,145 

2,041 

2,704 

2,381 

2,189 

2,739 

1,395 

2.227 
2,522 
2,812 

2.228 



Table II. Gross Weekly Inputs. Weeks 27 through 52, FY '74. 
(Continued) 
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Period 


analyzed - FY '74 


365 


days 


Shipments 




124,592 


shipments 


Volume 


of 


cargo 


10,070,724 


cubic feet 


Weight 


of 


cargo 


176,213,333 


pounds 


Number 


of 


PODs 


157 




Number 


of 


consignees 


3,356 





Table III. Data Summary for the Distribution Audit. 



SHIPMENTS 



Total shipments of 16 PODs 


115,732 shipments 


Total shipments of all PODs 


124,592 shipments 


Percentage of total shipments 


92.89% 


VOLUME 

Total volume of 16 PODs 


9,370,693 cubic feet 


Total volume of all PODs 


10,070,724 cubic feet 


Average volume/shipments (16 PODs) 


87 cubic feet 


Percentage of total volume (16 PODs) 


93.05% 


WEIGHT 

Total weight of 16 PODs 


160,575,726 pounds 


Total weight of all PODs 


176,213,333 pounds 


Average weight/shipments (16 PODs) 


1,387.47 pounds 


Percentage of total weights (16 PODs) 


91% 



Table IV. Tabular Data for 16 Major PODs , FY ’74. 
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POD 




SHIPMENTS 


VOLUME 


WEIGHT 


RA3 




20,125 


1,474,747 


28,284,704 


RGU 




6,041 


491,465 


11,486,615 


RG1 




1,850 


170,766 


3,774,348 


RJ1 




2,848 


170,635 


3,323,333 


RJ3 




2,210 


106,150 


1,802,072 


SA3 




16,837 


1,377,838 


19,920,384 


TA1 




2,502 


120,073 


2,132,222 


TA2 




6,793 


502,989 


7,269,693 


UB1 




12,339 


1,081,519 


19,311,360 


UC2 




6,593 


976,869 


19,325,856 


UD6 




8,691 


451,466 


9,062,837 


UL7 




3,183 


. 233,232 


4,445,639 


UM1 




6,619 


434,615 


7,081,348 


UM4 




3,769 


276,100 


3,740,452 


XE2 




14,008 


1,315,539 


16,679,435 


XJ1 




1,324 


186 ,690 


2,935,428 


Note : 


(1) 


Volume is in 


cubic feet 






(2) 


Weight is in 


pounds 





Table V. POD Data - 16 PODs , PY '74. 
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POD 




NUMBER OF 
RU SHIPMENTS 


PERCENTAGE OF 
RU SHIPMENTS 


RA3 




384 


1.90 


RGU 




156 


2.58 


RG1 




37 


2.00 


RJ1 




41 


1.43 


RJ3 




15 


0.68 


SA3 




183 


1.00 


TA1 




15 


0.59 


TA2 




86 


1.26 


UB1 




242 


1.96 


UC2 




327 


4.90 


UD6 




80 


0.92 


UL7 




43 


1.35 


UM1 




80 


0.92 


UM4 




20 


0.53 


XE2 




136 


0.97 


XJ1 




55 


4.15 




TOTAL 


1,900 




The percentage 
shipments for the 


of RU shipments 
16 major PODs is 


to the total number of 
1.64%. 



Table VI. Release Unit Shipments. 
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POD 

RA3 

RGU 

RG1 

RJ1 

RJ3 

SA3 

TA1 

TA2 

Note 



WEEK 
I NT 


WEEK 

START 


NO. d-STATISTIC 
OBS MEANING 


POD 


WEEK 

INT 


WEEK NO. d-STATISTIC 
START OBS MEANING 


1 


1 


52 


INCN 


UB1 


1 


1 


52 


NONE 


2 


1 


26 


INCN 




2 


1 


26 


NONE 


2 


2 


25 


POSI 




2 


2 


25 


NONE 


1 


1 


52 


POSI 


UC2 


1 


1 


52 


NONE 


2 


1 


26 


INCN 




2 


1 


26 


INCN 


2 


2 


25 


INCN 




2 


2 


25 


INCN 


1 


1 


52 


POSI 


UD6 


1 


1 


52 


INCN 


2 


1 


26 


POSI 




2 


1 


26 


NONE 


2 


2 


25 


POSI 




2 


2 


25 


NONE 


1 


1 


52 


POSI 


UL7 


1 


1 


52 


POSI 


2 


1 


26 


NONE 




2 


1 


26 


POSI 


2 


2 


25 


NONE 




2 


2 


25 


POSI 


1 


1 


52 


NONE 


UM1 


1 


1 


52 


POSI 


2 


1 


26 


NONE 




2 


1 


26 


NONE 


2 


2 


25 


NONE 




2 


2 


25 


POSI 


1 


1 


52 


POSI 


UM4 


1 


1 


52 


POSI 


2 


1 


26 


INCN 




2 


1 


26 


NONE 


2 


2 


25 


POSI 




2 


2 


25 


POSI 


1 


1 


52 


INCN 


XE2 


1 


1 


52 


POSI 


2 


1 


26 


POSI 




2 


1 


26 


POSI 


2 


2 


25 


INCN 




2 


2 


25 


POSI 


1 


1 


52 


INCN 


XJ1 


1 


1 


52 


POSI 


2 


1 


26 


INCN 




2 


1 


26 


POSI 


2 


2 


25 


INCN 




2 


2 


25 


INCN 


: (1) POSI 


' = the test ind: 


Icated 


POSITIVE serial 


correlation 



(2) I NGN = the test was INCONCLUSIVE 

(3) NONE = the test indicated no serial correlation 



Table VII. Serial Correlation Test Results. 
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MADs for 16 PODs. 
One-week time horizon 



Table VIII. 



POD 


TWMA 




OMA 


EWA 


"EWA-TA 


EWA 


-ARR 






t 




a 




a 




I 




RA3 


6 79 7 : 


* 3 


6797* 


.36 


7021 


. 11 


7186 - 


. 78 


7074 


RGU 


3910 


6 


3733 


.27 


3599* 


. 13 


3717 


.40 


3633 


RG1 


1754 


10 


1467* 


.60 


1709 


. 34 


1787 


. 34 


1611 


RJ1 


1793 


9 


1462* 


.35 


1734 


. 14 


1770 


.01 


1604 


RJ 3 


1130 


8 


1064 


.21 


1083 


.11 


1066 - 


.01 


1030* 


SA3 


7862 


3 


7862 


.04 


7458* 


.02 


7461 - 


.05 


7974 


TA1 


996 


5 


811 


.22 


935 


.10 


961 


. 13 


395 


TA2 


3654 


10 


3498 


.40 


3505 


.18 


3469* - 


.02 


3484 


UB1 


5321 


7 


5040* 


.24 


5112 


. 10 


5174 - 


.85 


5216 


UC2 


5236 


7 


4583* 


. 10 


4940 


.04 


4839 - 


.17 


4794 


UD6 


4205 


7 


3718 


• 

t'O 

oo 


3845 


.11 


3893 - 


.07 


3372* 


UL7 


1383 


5 


1296 


.20 


1290 


.12 


1358 


.38 


1287* 


UM1 


3637 


10 


3341 


.54 


3300 


. 10 


3311 - 


. 10 


3161* 


UM4 


2521 


8 


2366 


OO 

• 


2380 


.21 


2490 


.11 


2355* 


XE2 


7770 


10 


7391 


.47 


7615 


.26 


7787 


.07 


7130* 


XJ1 


2743 


2 


2564 


.43 


2434* 


.21 


2657 -1 


.00 


2475 


Note : 


(1) 


MAD 


= Me an 


Absolute 


Deviation 








(2) 


POD 


= Port 


of 


Debarkation 







(3) TWMA = Three Week Moving Average 

(4) OMA = Optimal Moving Average 

(5) t = time interval (in weeks) for OMA 

(6) EWA = Exponentially Weighted Average 

(7) a = smoothing constant for EWA and EWA-TA 

(8) EWA-TA = Exponentially Weighted Average with 

Trend Adjustment 

(9) EWA-ARR = Exponentially Weighted Average with 

Adaptive Response Rate 

(10) y = smoothing constant for EWA-ARR 

(11) MADs are expressed in cubic feet 

*Indicates smallest MAD for POD 



Table VIII. f-IADs for 16 PODs (Continued). 

Two-week time horizon 



POD 


TWMA 


t 


OMA 


a ; 


EWA 


a 


EWA-TA 


] 

I ' 


EWA-ARR 


RA3 


6204 


3 


6204 


. 32 


5955* 


. 17 


6293 - 


.89 


6095 


RGU 


3621 


3 


3621 


.09 


32.84* 


.22 


3489 


.45 


3339 


RG1 


1601 


10 


1362 


. 77 


1572 


. 37 


1690 


.23 


1311* 


RJ1 


1350 


8 


1194* 


. 37 


1420 


. 17 


1449 


.01 


1292 


RJ3 


850 


3 


850 


.28 


804 


. 11 


801 


.01 


750* 


SA3 


6639 


3 


6639 


.02 


6250 


.01 


6247* - 


.05 


6927 


TA1 


798 


6 


595* 


.22 


701 


. 12 


729 


.11 


694 


TA2 


3194 


10 


3056 


. 39 


3072 


. 17 


3125 - 


. 03 


2991* 


UB1 


4655 


7 


4250 


. 15 


4246 


. 06 


4179* - 


.01 


4417 


UC2 


4001 


7 


3150* 


. 15 


3399 


. 06 


3278 - 


. 10 


3423 


UD6 


3688 


7 


3135 


.27 


3235 


. 10 


3268 - 


.08 


2938* 


UL7 


1053 


4 


1005* 


.28 


1013 


. 17 


1082 - 


.44 


1033 


UM1 


3146 


10 


2641 


. 19 


2816 


. 09 


2893 


.01 


2615 


UM4 


2343 


9 


2153 - 


.02 


1878 - 


.01 


1847* 


. 10 


2177 


XE2 


6984 


3 


6984 


.65 


6697 


. 31 


6941 - 


.01 


6269* 


XJ1 


2771 


10 


2327* 


.23 


2435 


. 10 


2592 - 


.75 


2452 



indicates smallest MAD for POD 



Table VIII 



MADs for 16 PODs (Continued) . 



Three -week time horizon 



POD 


TWMA 


t 


OMA 


a 


EWA 


EWA-TA 

a 


EWA-ARR 

Y 


RA3 


6620 


8 


6197 


.26 


5910 


.10 


6069 


- .74 


6351 


RGU 


3477 


5 


3410 


. 13 


3120* 


.04 


3382 


- .47 


3214 


RG1 


1652 


10 


1363 


. 75 


1498 


. 37 


1529 


. 25 


1245* 


RJ1 


1188 


7 


1036* 


.43 


1253 


.25 


1296 


.01 


1086 


RJ 3 


863 


10 


824 


. 31 


789 


.09 


784 


. 01 


713* 


SA3 


6451 


7 


6388 


. 34 


5848 


.01 


5238* 


-1.00 


6234 


TA1 


721 


5 


545* 


.25 


626 


.13 


653 


. 14 


561 


TA2 


3149 


10 


2660* 


.31 


3030 


. 17 


2990 


- .03 


2789 


UB1 


3627 


7 


3252 


. 18 


3225 


.07 


3187* 


1 

O 
1 — 1 


3471 


UC2 


3390 


10 


2548* 


. 19 


2830 


. 06 


2781 


o 

o 

T— 1 
1 


2652 


UD6 


3326 


10 


2981 


.29 


2988 


. 11 


3147 


- .07 


2410* 


UL7 


972 


5 


937* 


.29 


946 


. 16 


995 


- .46 


978 


UM1 


3215 


10 


2355 


. 19 


2675 


.09 


2681 


- .12 


2331* 


UM4 


2327 


9 


1970 - 


.01 


1709* 


.01 


2051 


- .01 


2127 


XE2 


6831 


10 


6717 


. 55 


6778 


.30 


6960 


- .01 


5971* 


XJ1 


2787 


10 


2049* 


.17 


2308 


.09 


2422 


. 75 


2386 



^Indicates smallest MAD for POD 
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Table IX. MDSs for 16 PODs 
One-week time horizon 



POD 


TWMA 


t 


OMA 


a 


EWA 


a 


EWA-TA 


I 


EWA-ARR 


RA3 


92,564 


7 


90,250* 


. 30 


93,854 


. 06 


98,367 


- .82 


98,528 


RGU 


27,737 


4 


27,546 


.32 


24,977* 


. 14 


26,690 


- .59 


25,060 


RG1 


10,349 


10 


5,769* 


.60 


10,736 


. 30 


10,117 


. 25 


8,805 


RJ1 


5,741 


8 


3,691* 


. 33 


5,337 


.16 


5,511 


.02 


5,474 


RJ3 


2,825 


7 


2,591 


.12 


2,455* 


.05 


2,481 


- .55 


2,455 


SA3 


96,136* 


3 


96,136* 


. 38 


98,728 


. 04 


304,177 


1 

I— 1 

O 

O 


104,207 


TA1 


1.659 


6 


1,184* 


.21 


1,545 


.12 


1,611 


. 11 


1,614 


TA2 


22,695 


10 


19,113* 


.28 


20,456 


. 12 


21,253 


- .03 


20,549 


UB1 


51,088 


7 


41,265* 


.14 


42,697 


. 06 


43,547 


- .78 


44,880 


UC2 


41,418 


9 


30,695* 


. 14 


35 ,7 75 


.05 


34,811 


O 

O 

t-H 


35,382 


UD6 


34,556 


7 


30,086 


.25 


29,872 


. 13 


31,267 


- .07 


27,730* 


UL7 


3,453 


5 


3,195 


.28 


3,183 


. 16 


3,267 


- .44 


3,176* 


UM1 


25,750 


10 


23,983 


.24 


22,108 


. 11 


23,010 


- .14 


21,463* 


UM4 


11,870 


8 


9,145* 


.28 


9,936 


.07 


10,610 


- .86 


10,500 


XE2 


102,515 


10 


84,346* 


. 52 


103,848 


. 25 


106,797 


.05 


90,932 


XJ1 


15,387 


2 


13,530 


.48 


12,845 


. 16 


14,115 


- .13 


12,490* 



Note: (1) MDS = Mean Deviation Squared 

(2) POD = Port of Debarkation 

(3) TWMA = Three Week Moving Average 

(4) OMA = Optimal Moving Average 

(5) t = time interval (in weeks) for OMA 

(6) EWA = Exponentially Weighted Average 

(7) a = smoothing constant for EWA and EWA-TA 

(8) EWA-TA = Exponentially Weighted Average with 

Trend Adjustment 

(9) EWA-ARR = Exponentially Weighted Average with 

Adaptive Response Rate 

(10) y = smoothing constant for EWA-ARR 

(11) MDSs are in 1,000's and expressed in (cubic feet) 2 
^Indicates smallest MDS for POD 
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Table IX. MDSs for 16 PODs 
(continued) 

Two-week time horizon 



POD 


TWMA 


t 


OMA 


a 


EWA 


a 


EWA-TA 


X 


EWA-ARR 


RA3 


70,264 


7 


69,626 


.18 


62,116* 


.05 


64,218 


- .82 


69,925 


RGU 


21,534 


5 


20,533 


.27 


18,672* 


. 12 


20,411 


- .53 


18,829 


RG1 


9,064 


10 


5,293* 


.69 


8,385 


.36 


7,121 


.22 


5,588 


RJ1 


3,732 


8 


2,413* 


. 32 


4,019 


. 16 


4,165 


.01 


3,965 


RJ 3 


1,659 


10 


1,388 


. 13 


1,327* 


. 06 


1,360 


- .55 


1,334 


SA3 


76,062 


3 


76,062 


.07 


73,856* 


.03 


73,992 


l 

I- 1 

O 

o 


80,384 


TA1 


1,060 


6 


749* 


.22 


953 


. 13 


1,018 


. 11 


948 


TA2 


16,801 


10 


12,717* 


. 26 


14,376 


. 11 


15,107 


- .03 


13,101 


UB1 


30.872 


7 


24.900* 


. 14 


25,001 


.06 


25,749 


- .01 


27,513 


UC2 


25,069 


10 


14,927* 


. 15 


18,843 


. 06 


17,843 


1 

h- 1 

o 

o 


17,496 


UD6 


24,416 


10 


20,683 


. 23 


21,071 


. 12 


22,598 


- .07 


17,668* 


UL7 


2,084 


4 


1,880* 


.30 


1,952 


. 17 


2,021 


- .47 


1,977 


UM1 


19,208 


10 


15,958 


.20 


16,063 


. 10 


16,948 


- .01 


14,948* 


UM4 


10,062 


9 


6,781* 


.15 


7,909 


.01 


8,057 


- .01 


8,216 


XE2 


81,399 


2 


80,557 


. 54 


79,063 


. 26 


81,152 


- .01 


60,822* 


XJ1 


13,673 


10 


10,908* 


. 31 


11,087 


.13 


12,110 


- .14 


11,067 



^Indicates smallest MDS for POD 



Table IX. MDSs for 16 PODs 
(continued) 

Three-week time horizon 



POD 


TWMA 


t 


OMA 


a 


EWA 


a 


EWA-TA 


I 


EWA-ARR 


RA3 


69,762 


8 


64,618 


.10 


54,579* 


.04 


55,583 


- .78 


67,596 


RGU 


.19,006 


4 


18,391 


. 26 


16,418* 


. 11 


18,245 


- .51 


16,570 


RG1 


10,187 


10 


5,461* 


. 75 


8,571 


. 38 


6,724 


.24 


5,822 


RJ1 


2,737 


10 


1,871* 


.33 


3,234 


. 18 


3,354 


1 

O 

h-* 


2,968 


RJ3 


1,347 


9 


1,035 


.14 


967* 


. 06 


1,002 


- .55 


969 


SA3 


74,786 


3 


74,786 


.05 


61,446 


.02 


61,394* 


-1.00 


73,494 


TA1 


806 


10 


663 


. 26 


712 


.14 


77 9 


.11 


635* 


TA2 


15,444 


10 


9,783* 


.24 


12,773 


. 10 


13,408 


- .03 


10,545 


UB1 


18,965 


4 


16,064 


. 17 


15,866* 


. 07 


16,757 


- .01 


17,605 


UC2 


19,213 


10 


10,170* 


. 16 


13,467 


. 06 


12,539 


1 

O 

O 


11,723 


UD6 


20,080 


10 


15,516 


.24 


16,499 


. 13 


18,160 


- .07 


12,549* 


UL7 


1,633 


4 


1,551* 


.33 


1,570 


.18 


1,638 


.49 


1,621 


UM1 


17,862 


10 


12,662 


. 19 


13,577 


.10 


14,584 


- .12 


11,984* 


UM4 


8,787 


8 


5,197* 


.01 


6,067 


. 01 


6,217 


- .69 


7,0 84 


XE2 


81,188 


10 


75,150 


. 53 


79,907 


. 26 


80,197 


- .01 


54,919* 


XJ1 


13,778 


10 


8,106* 


.23 


10,803 


. 12 


11,790 


- .94 


10,820 



^Indicates smallest MDS for POD 
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WEEKLY 

TIME 

HORIZON 


TWMA 


t 


OMA 


a 


EWA 


a 


EWA-TA 


Y 


ARR 


one 


545,749 


7 


523,120 


.37 


535,979 


. 17 


566,904 


- . 14 


543,745 


two 


406,966 


10 


391,469 


.36 


388,940 


.17 


421,432 


- .05 


393,060 


three 


375,569 


10 


341,400 


. 36 


350,451 


. 17 


384,810 


- .05 


346,852 



Note : 

(1) MDS = Mean Deviation Squared 

(2) OMA = Optimal Moving Average (a specific optimal interval 
is applied to all POD data) 

(3) EWA = Exponentially Weighted Average (a specific optimal 
a-level is applied to all POD data 

(4) EWA-TA = Exponentially Weighted Average with Trend Ad- 
justment (a specific a-level is applied to all POD data) 

(5) ARR = Adaptive Response Rate (a specific optimal y-levd 
is applied to all data 

(6) t = optimal time interval for OMA 

(7) a = smoothing constant for EWA and EWA-TA 

(8) y = smoothing constant for ARR 

(9) MDSs are expressed in (cubic feet) 2 



Table XI. Aggregate MDSs for 16 PODs . 



In this work the Student t-distribution was used to 



resolve whether the mean of the MAD differences between two 
forecasting methods was equal to zero. 

The hypotheses were: 

H q : y = 0 , H x : y f 0. 

The t-test utilizes the following 

t = U - u) /n 
s 

where X is the sample mean, y = 0, n=number of observations 
and s=the sample standard deviation. 

Calculated values of t are compared with critical values 
of t. If 

t (calculated) < t (critical) 

H q is accepted 

If t (calculated) > t (critical) 

H q is not accepted. 

If H q is not accepted, the two methods are different and 
the method that generates the lowest mean absolute deviation 
is chosen as the preferred. 

The test only utilized deviations generated by a one- 
week time horizon. 



Table XII. t-Test Procedure. 



COMPUTER OUTPUT ONE 



I 

TABLE Of STATISTICAL DATA FOR SFLfcCTcD PODS USING WEEKLY TOTAL VOLUME (CU. FT.) DATA 



POD 


INIEKVAL 


START 


INTERCEPT 


SLOPE 


ST u, E KR • 


R-S0UAR6 


D-STAT 


OBS. 


ME ANIf 




I 


1 


0. 


0 . A 1 Ao 7E OO 


l 75. 1 


0.00131 


1 .60802 


52 


MU ME 




2 


1 


125. 


0.1781 9 E 01 


2A7.5 


0 • 003 1 5 


1 .3A05 8 


26 


IN CM 




2 


2 


13. 


0 • 82 70 If 01 


239.2 


0.06791 


1 . 30oQ 6 


25 


NONE 


AC 2 


I 


1 


2. 


0.99676E-01 


28. A 


0.00289 


2.07288 


52 


NON r 




2 


1 


5. 


0.35385E 0 J 


AO. 2 


O.OOA69 


2. IA859 


26 


NONE 




2 


2 


7. 


0 • 3 1 1 7gE 00 


h0.3 


0.00363 


2.1652A 


25 


NONE 


A FA 


I 


I 


56. 


-0.11562E 01 


133.9' 


0.01717 


1 . 99 7 AO 


52 


NONE 




2 


1 


115. 


-3. A8520E 01 


167.7 


0. 03WA 


l . 1 163A 


26 


P0S1 




2 


2 


AO. 


-J.AdAl l E -0 1 


183.2 


0.00223 


0.98535 


? 5 


POSI 


BA 1 


1 


1 


973. 


-0.1A17AE 01 


1588.2 


0.00019 


1 . 79631 


52 


NONE 




2 


1 


1936. 


-O.A8l23t 01 


2A32.2 


0.0002** 


1 .751A5 


26 


NONE 




2 


2 


1901. 


-J.66281E 01 


23oO.B 


0.05078 


l .8 20 76 


25 


NONE 


BA 3 


I 


I 


3. 


-0.576c BE - 0 1 


12.6 


0.00A90 


2 . 0 A 9 2 7 


52 


NONE 




2 


1 


6 . 


-0.21 535 E 00 


l 7.9 


0.00071 


2 .0 y d *♦ 3 


2T> 


NUNF 




2 


2 


220. 


-0. 1281 Ufc 02 


277.3 


0.11508 


0. 1 1955 


25 


POSi 


BAA 


I 


I 


1 B A . 


-0.1136 IE 01 


393. 1 


0.00155 


l . 69 A09 


52 


NONE 




2 


1 


3/3. 


-O.A7521E 01 


5 2 A . 6 
621.2 


O.OOA98 


1.68739 


26 


NONE 




2 


2 


•A A 1 . 


- 0.13197b 32 


0.026/7 


2 . 1 A 356 


25 


NONE 


CA2 


1 


1 


-9. 


0.5A517E 00 


ao. e 


0 • 0 A 0 2 A 


2.12489 


52 


NONE 




2 


1 


-17. 


0 . 2 1 3 7 IE 01 


57.1 


0 . 0 7 8 A 0 


2 . 2 58 AO 


26 


PUNE 




2 


2 


AA. 


- 0 . 1 3 1 8 * c 01 


101.7 


0 . 0 1 0 1 A 


1 . 7A0 70 


25 


NONE 


CA3 


1 


1 


1. 


-0 . 2 1 5 l 5E-0 1 


3.3 


0. 039o0 


2 .Os j 02 


52 


NONE 




2 


1 


2. 


-0.90 2366- 01 


3.8 


0.02151 


2.1261 1 


26 


NuNE 




2 


2 


3. 


-J.AA162E 00 


9. 1 


0. 12650 


1.55352 


25 


NONE 


C65 


1 


1 


10. 


-0.375o5E-01 


50.6 


0.00013 


2.07061 


5 ? 


N jNE 




2 


1 


21. 


- 0 . 1 005 IE 0 j 


/i. a 


0.00035 


2. !Ao62 


23 


NONE 




2 


2 


19. 


-0.20310E-0L 


71 . 7 


0 • OOOuO 


2 • 1 5 a? 1 


25 


<koNE 


CL I 


1 


1 


6. 


0.?3h 36t-ul 


38.3 


0 .00009 


2.0526A 


5 2 


NONE 




2 


I 


13. 


0.35500E-01 


5 A . A 


0.00003 


2.1 1631 


26 


NONE 




2 


2 


1A. 


-0 . 1 1 74 6 E 00 


5*t. A 


0.00020 


2.113/5 


25 


NONE 


CE 2 


1 


1 


-0. 


0.A9051E-02 


0.7 


0.01151 


2.06302 


52 


NONE 




2 


I 


-0. 


0.18800E-01 


1 • u 


0.02151 


2.12611 


26 


NONE 




2 


2 


2. 


-0.99167E-01 


2.5 


0.06/67 


1 • 0 3B A 8 


25 


POSI 


CHI 


1 


1 


-1. 


0.66166E-01 


6.9 


0.02091 


2.08287 


52 


NONE 


2 


I r 


-2. 


0.256AIE 00 


9.8 


O.OAOOO 


2.16737 


26 


NONE 




2 


2 


- 1. 


0.27A3AE 00 


9.8 


0 . 0Ai>9 3 


2.20913 


25 


NONE 


CJ1 


I 


1 


0. 


0.57629E-02 


2. 1 


0.00176 


2.0A282 


52 


NONE 


2 


1 


0. 


0.236Alt-0l 


3. 0 


0 • 004 At 


2 .08936 


2 6 


NONE 




2 


2 


0. 


0 . 35371E-01 


3.0 


0 . 0 0 0 A 5 


2.0723A 


25 


NONE 


CKA 


I 


1 


27 5. 


-0.A27A7E 01 


39 A. 5 


0.0 25 7 7 


l .A155A 


52 


POSI 


2 


1 


555. 


-U.17ls9= 02 


572.3 


0.05194 


1 .60387 


25 


NONE 




2 


2 


3 A 0 . 


-U.16553E 01 


713.3 


0. 01205 


2.0561 2 


25 


NONE 


LK1 


l 


1 


7 1. 


-0.11 2 AGE 01 


129. A 


0.0173/ 


1 • i 40 A 1 


2 


POSI 


2 


l 


1 A . 


-J.A555oc ul 


18 5.9 


0.0338V 


l . 1 l 1 79 


^6 


POSI 




2 


2 


1A5. 


-0.A7130E 01 


237. 3 


0.02^7 


1 .5A2d5 


25 


NONE 


CK 2 


1 


1 


135. 


-0 • 30 8 7 JE Jl 


278. 8 


0.‘02 /93 


2. 05219 


d2 


NONE 


2 


1 


268. 


- 0. 1 20 1 6t 02 


385.5 


0 .05A81 


2 . 18693 


26 


NuNE 




2 


2 


27s. 


-O.12970E 02 


AO 3 • A 


0 .05927 


2.26802 


25 


NONE 


EOl 


1 


1 


0. 


-0.73423E-0? 


0.5 


0.040?A 


2.1 ? A p > 7 


52 


NONE 


2 


1 


1. 


-0.2871 BE-Ol 


0.8 


0. O/oAO 


2 . 2 5 8 1 0 


2 6 


N <j N t 




2 


2 


38. 


-0 .2251 7 L 01 


A 5 • 0 


0. 11A11 


0.06722 


25 


POSI 


ED 3 


1 


l 


0. 


-0.3158 9 F-02 


0.3 


J. 02979 


2.10197 


57 


NONE 


2 


I 


0. 


-0. 12991E-01 


0. A 


0.06** 18 


2.22332 


26 


NONE 




2 


2 


5 « 


-0.319/&6 UJ 


6 • C 


0. 11932 


0 . 1 l Ac 3 


25 


PUSi 


EOA 


1 


1 


25. 


-0.36933E 00 


32.1 


0.03002 


1.83209 


52 


NONc 


2 


1 


50. 


-0.1A988E 01 


A6 . 6 


0.05923 


1 .6A072 


26 


NONE 




2 


2 


A6. 


-0.12378E 01 


50.1 


0.03582 


1.A8792 


25 


NUNE 



101 



COMPUTER OUTPUT TWO 



► ***#*< 



**ftft*<}**#» ft * a * < 






Csigr- 

rA A <A g •— g ajN-r»r» 
00 f-0)-4 rr^rocoosorn 

0 — «rgAgN» -o — iaa -o 

1 OOMNifiOO'OO'l 
Ag xiNeciNin*;^ o 
n* rt> rsj rg eg eg rg rg eg rg 



rg «0 OOOfNJrg t)'0>0 
As/'r- oogr-o t <7* 
— «Njg(T).Ag *Q'Y) T) 0 
mg g g in orgiron 
OAg-*rg,A oaos— < 
g * 0) 

O' co rg eg rg eg o — • O' O 
rgrgrgrgfgrgogrg^^o 



m — • g ro 
r^O'O'O 
>fON^ 

rocuroo 
o — *a — « 

egggrg 
A rOrO "A 



grorg ao o 

jsinoop'jj 
nini^inrgo 
go AA>otA 
gegg ojg 
(ArO rAcOrorA 



O'JrgrD 
— *'\j-org 
O OO-* 
ao reoeg 
ao 
rgog o 
r-o 
OAgg 



>00 0 r rgoi 
go r o rr rn 

g o gg nj> 
■or>* ur^r-n 
g r 03 <A g— < 

'fiftf’tN O CO 
— <000 
rA*0‘ArAC0r0 



O g— «ON~gr-oggA A A g 0r003 g g coco 

A rO rA OiONN O ON- r"-OA — 

■“ * — OJ003 OAr- oN-Ag 



an- mg 03mgcgx>r , ~ 
■onoo3 r gog ggg 
N~ >r m rA egm eg ro rA g 
— grAmrAmmrAmmrM hm«mh - » ih. < ■ i>h 



03cg03 f '-majimAmr- 

mg ,a J'-joo'odo om 
MnfNJfOO'iTuTlO OO 
03 ausus vr g g g g g 



O — «rg-ogAOr^aogO 
>■ —* 
< 



^rg(AgA,or»aogo — «rgmgA 03 r'-aogo — 'rgrAgAO^r— g)go 



o 

o 



»A 

3 

a 

♦ ftftft*ttft0*ft 
*ftftft>«ft«ftft 



ftttftftftttftfcftft 

ftftftftftftft*ft# 



ft ft ft ft ft ft ft ft ft ft 
ftftftftftftft*ftft 



rg 
O 
. J> 

ft ft ft ft ft ft ft ft ft ft 
ftftftftftftftftftft 



O 0Ng‘33(T>f07 0'^ 

gcogr^mg >ton g 

03 03^U3gc0ggga0 
OrAgcGAi'O— *«3g J* 

j-g-gnNNr »g 03 

ng-^-OMco-n^NN 

cogog-^rnggr-A 



OA— 'Org.gg-rgco 
g g a com g ro A g ^ 
rgu'i— <rg-oOAr"tng 
g o g <\j gg 0000 

fOO^rgiga)^f\i«on 

oj — « 03 A»\j— • — * rA<\j* »» 



O orgro g<M 00 0<g 
wOrnr^ogogN. g gmg 
onons doo^mti 
o 

gr»mcgir\g ogusog 
O'OUSiO ^(OOOON 
rgm ^HfM^oooo 
gcofAmmmmmrom 



0003-0 -roajrgg 
— * 03g -OrgegcorgmAj 
OAO — *r- g 03 a -org 
eg OvTM'liornoO'O 
>ou— nnor*inusNg 
— • a a g a o g — < Oca 
-tmrgooo-^fMocug 

fOOO-'Hisj^uu'io 



O fO^gNtn-HivgNN 
< 'OAACOO g g> 

i: mAr«AN« cor^ 00 g 



w-« 03 rg-H^gogr-co NrogSrsusogiNON 



03— «g g— «rgrgcoA N- 
rg^gaocoaocococDco 



O' rg o eg — » g- — 1 a O' A 
03 g iNJ g CO 00 r- N- N- N> 
g •ggmmmmmmm 



rVilOONf\J^rjJ^^_4 

rAN-r-cQrAgr>-omg 
gr^r~ g— «o an- 03 m 
e-r-r-N-mJur-r'-r-N. 



CD ~-«rvjmgA<ON*mgo — •rgmgA'ON-aoOO 
> — ♦ -< 
< 



-Hrgrogmg)N-oo< 7'0 — ♦rgroggsgjf^aoO'O 



CD 

UJ 

o 



> 

-t 



O 



o 

£ 



. #ftftft*ftftftft* 



ftftftftftftft*ft# 

ftftftftftftftftft# 



o 



ft ft « ft ft ft ft ft ft * 
ftftftftftft aft « ft 



ftftftftftftw-ftft# 

ftftUftSftft-cftft 



OOO geo OO (M 03 o 
Nirioorgoo<;C'g3g 
snoDODrorga3(j'JDaoo 
O'ONNgjgogggs'o 
rsgooooDgoDcoo3fv*g 
-Hsrs^'X)srsg'gso3ror s ' 
m<7'N-r~r~r-fT'^-«f\jrg 
fOrgrgrgrg >g . gn.-nro 



0000 ■O'gorgrgg 03 
rg rg -h rn m OJ r~ 03 OD 
Nmotno^ONOg- 
nooMnNjsjsooiri 
■ocrrooaorgrg U'J'O 
rO UO — < Ocn — *0 U 3 X 03 

noo og n(Mgina3i3s 
rg— <jso— • og — « 1^- i^-4 



O' g no -4N- co g g ro 
jooja'r^romgc— *T) o 
*0 co x) O' — «rg oc'N-g 
OONgrslH^oO® 
CD cO 03 O r^J r — < 3 -A 
OD com Oig^OgON 

g«nurs^n oonnn® 



g3— iOr^oDin— *>00 
ONiS -o OCTN-riajro 
■ggotnooosgoo 
— I OOO'OUNJ'® gcr* 
or-r~nc7 u-o-*gm</> 
•ooor^vrtg-^r'.-g 
rorg— *— <OoCTs7'C> f 7' 



2 

U 



rSrnog-Oro— <rgO— * 
0)00— <<gcoro<g o3sfs^ 
1 sn rg^U'aoN* 0000 
ggrsArmvnggg 



O t fg^NLDcOcogN 
CTiTi 03 — <N-Orgao — <«-h 
oenco— «g ,, -'n— <^0,0 
'7'aDN-aDOOa3DOCOcoa3 



co g.'O^orgN-uoaDN- 
iA0)ga0uD(0N-0D03O 
Ojinr-r^'oa' J'Ocr x> 
rgrgrgrgrg^rg nrg rg 



O— f-^OirDOm 03 og^ 
m rg r g o 0) co n- oo eg o 
r-insnoDggr*3,ng g 
r g rg rg 1 g rg rg rg cnj rg rg 



O — 4 <M" 3 gA 03 N- ®NO 

> -« 



-•rgrogsT. Or^TDO'O — irg^gg^ o)r-coo s O — ♦rg.t^ gg\ or-u>C'0 



CD 

2 



O 

O 

a 



3 

o 

ac 

ftft»ft»ftftftftft 
ftftftft ftftftftftft 



ftftftftftftftftftft 
ft ft ft ft ft ft » r. « » 



•ft ft ft ft ft ft ft ft ft ft 
ftftKftftftftftftft 



3 



ftftftft ft ftftftft ft 
ftftftftttftftftftft 



O 

o 



JC 

o 



2 

o 



rgrgo3ggrgo)03Ajg 
fOr*cr>ajaj— 4*^.o^u^o3 
co 03 gNi-ftNNg .no 
rgOgoO®NysrriN x> 
mg oN-gmgrgog 
r— ojosr^N-— «cgr^moj 
Oenrgg^vT'Orgrog 
rocftcra'oa'cra'cr'a' 



NvogNuonu'^^'O 
cjugsgesroiS-NioDino 
O' 03 rg — *n- n- — • g in 03 
com— * qnnosh'Du'S 
^ggO®®r\i7gu 
mo)N'gorgN.'UNO‘ 
n. omgggfDfftrftri 



f^NOW 

tn'ooo 
COg X; 4S 
g<\n3— < 
03 00 rN 

m 03 — < g 
n - 03 sng 



O r 0) OOrg 

"M r — • 'U-t 
■OlNrjNg-ft 
f'- vu n g^N 
ocNo;inrgo3 
OD-Nrg 3 0)0) 
gos-tiMftig 

gggggg 



0 O uo a3 1 M O> g vX> 03 O 
rg— <0 >r J' O-uO 03 co 

osgno^N oogrg 
rorgo o— • O— •— iN,-n 
g r 1 u s cn rg ^ co 

g oj n- r- g n» 0 3) g 

g03 0ftroiA03 03ggrr> 

1 g eg --g rg rg nj f\i eg rg 



q g^r-rorgg ^ooftrg co oDrnrgrAOft.niArgoft 

<J yMTlO'Ng-ogNN® (jiHU'NN'" *NC0 03— • 

71 trs 03 N» 03— g g gr«-iAgir> gg g'A 
aor- ojo3N-n-03N-n-n~ — 



CD 0-*0>OgOsA-isO 
iA jdosj -\j rgo g co ro 
o fA rg Njrrj o — • ro g 
r— »OiAiAsAiAiAsAsAiA 



rV 7 N cO rO N«- ro A C 23 — * 
OgrA0303 03iAgr^O 
N-fA03 ro 03 N- O-rrftH 
(A r3 fA ro rO ro fA ro "ft fr> 



> 
ui 
O 

U- o — «rgrrig sA>ONa3go — «rgrr»g A03N-cogo — *rgrA g a 03 n- ao g O -orgi<)g®>ON®^0 



a» 030000 jrggrgrgro 
o— *orgr*grgN-g g 
33 03 0 ggr^g co -Njrs 3 
CJO OOSmnrgf^g-Hr. 

?; od 03 g -HrAm— <cncn o 
On»r^r^.g <S — • coco r ) 

D rAoggr-CDrgrO'Arg 
J r~ ODt-TSiAAsAsAcAAsA 



L) 

I — •••••••••• 

O -ON-rg 03 N- g gr^rg— • 
LU<jgN~— <— <— 'rg'jcog— < 
h-ft. Ajor-gr-— »rgiA -o A 
*i oi'sogcj— »gu'v'g 
O N- 0)o3A 0) mJ AUS A A 



CJ m 

CJ <t 



co 

3 



X. 

3 



oo— *cg*A g A 'ONu/vro 
o> — « 

<J< 



102 



COMPUTER OUTPUT THREE -ONE 
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Smoothing Constant (alpha) 



COMPUTER OUTPUT THREE -TWO 
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COMPUTER OUTPUT THREE -THREE 
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COMPUTER OUTPUT FOUR-ONE 
Mean Absolute Deviation (in cubic 
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Smoothing Constant (alpha) 



COMPUTER OUTPUT FOUR-TWO 



« 



u 



at flo©'OD' 4 , >onO'Tin'r 'j , ooi^o-or v -i/ ir, nny' -o/grg m^^r- u*a>mo''gj cuao»-g-Htr o mm O' oo'Oia - 0 ^ 0 '^^ ou<rnnua''To o ond o^o 

4/1 m or-omr- omiA 'vjrg'U 4 rgiA^g^giAOOrg 4 <\jr- r 4}0~« >oo'(\j~«r-rgr , ~ rgmr-o> mmo* 4 r-r- mmcoo'<— • 4 r-- o 40 ' omrg j— imo* 

> ■nj'N mmo>o coN'D^Of'*TJ( , n<ofMin m y mo 'O ucgunmggrgrgm 4 y gjiA<— «m'-*ir\m 4 U'OAgjr-'Yjo' 00 ' O — « 

uj gJO'm xmo' 4m ~«rgmoiA .or— msotnorg ooujo mor-m r -iNin^M^^fT'iDo^^sr Din^ oof^'ON^Nir\iANH't)ininr»oiANOO'^ t'-O'f 



O' 'f '*«'or~'-O r — co^4com.x>4'-<4r— yyy.o-.no'''-J 4 u~«iAO'AO'-o~'JO'ymoayr— A 4 mrgc\JAjmmm 4 iAr-aoocg 4 r- O'miAcom n O' mr—^iAO' 4 
^ 'X’ 4 ^ooia mur— <r»rgo -o r-mro'-^ooo ■oiAmrgoo'r— o 4 mrg-go'aor— •oin 4 mcg^oo'cor- 4 jLAg \4 mm-g~<O0'0'x/r-r>- oyiA 4 t 
mmm <\jmr\j — «^^-hOOooo* 0 ' o* O' 'aox>aox>''oaor-r-r-r— y yyoy-o.oiAiAiA'A<AiA'AiAiA , A4 44’444444444mmmmmmmmmm 



> 

UJ y04cr4CT'U^4mAjr\jrgm4 y r- J'l'jincorgr-^rg 4 r- 0'f\J4comiAao*- • Acomr—cgr— rgujmor—mo' y mor- 4~*4'coTor—OA4r\JOO''— 4 mo O' O' 
lj crymy 4 *~«®mmo'ymor- 4 ~«oymocDAm^: 7 'r»m 4 moo'r— iA 4 moo'*—y 4 mmou'r-y 4 mm~« 0 'cor**iA 4 mm«"-*oo'cor-v^mig'-«(?'ao 

< 'J' v?- 4 mmmmmm~<— « ~*~*00 00 ' 0 ' 0 ' 0 'jocococor-r— r— r—r-r— .y -O'O'OO'U^^^'AinininNrr ^ 44 -4" 4mmmmmmmmmmmm rnmmmrn-*-* 

J. mmm mmm mmmmmmmmmm-«^~«^''-'--« -g-g ~*-g -g~* m 

< 

I OO'acr— '■OiA4mm-*00'aor-«OiA4mcg^OO'aor- yLA4mm ~*0 0' JOf— y iA4 mm^O oaor- OcA4mm^OO'X>r-yiA4mm”"<00'aor—yiA4mm 
o. oo' 00 '( 7 'V' 0 > C'C 7 'V'( 7 'cD'ygoajcDcoaoajcocDr-r s -r-r-r— r-h-r«-r— r-yyyyyyyyy 'OininirunLmrunmioin^'rvj- 44444 44 mmmmmmmrn 

< -lOOOOOOUOOOOOOOOOOOOOOCOOOOOOOO OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

o m 

o -> 

O. a 



O CO 4 CO -O'OOOO^OMOaOO-J'COOOO Uf\JCDO(N' 0 'D 00 > 0 C 3 C 0 T 3 'J , vrO' 0 a 0 O''JN' 0 ' 4 1 -^'/MAmv}-NNn^>J’> 0 N(\OOO-^O' 0 SC 0 ' 0 -«' 0 r^NN>t 
to y yoiA~«iAy cr'm 4 iAr-rg'UJOO~gy 4 y o' co riAr-oaoo'yoy 4 4 mo'OOtAr'-r-om 4 ~*mmm-*iAr-mma‘>cur-o'maooautAj' 4 C'mujr- 4 m 
> inmym04^0ym— (O a o — < 4 <x>«-^miAO' 4y^g^grgmmr— cocAm-A 4 4 mao j'Or'iAO’T'OD f 4r- C'O' itin mmiAmiA ror-in— <r-o t r- om y O' 

UJ m AO yO 0 'rgy'’ 0 C 0 mOOCr'r- '^O'fOr^sf^mNsr >OM'H'Mir>7'(<1vO<t)DOOHa' 4 rg 4) O' 4 m m r— 4 4 O —•COr-O'mOOm y y~ «ommaDm 4 O 

o oo'A«o 4 yiAco 'OU'n (mt> 400 yr-^cj' tTmocr^y mm Aa*iAm 4 oo o 4 4iAr— mr— mmoiA'X) mr— 4 ^rvjtnoo ym^oomy ~»r- la 4 o O' 4 
a. »•*<' >r- — » y ^ r— y r-com aqo .-* 4 mm y o*ao' -r <4 <r o* la y rg m 4 1 _? r- iior~ j-~-«yt.nmoyy 4 (■\r-J‘yr-A 4 mm'- , ~* 

'f'NO^N y 4 mmoo'OJf— y la 4 mm^oo'aor—r- y la4 4 mmm-*'-*O 0 ' 0 ' coaor— r-r— y y u\ala 4 ^vj-tnintniNMNfMH^H^^ouuoyooo 
iaaia 444 -J' ^ 'J - 4 m m m mm m mA m mm 



> 

uj co-^colA 4mmm 4LAaom'y 40f'-LAmm'-^— «'-gm<Ay^y — < OmcuiAmOaoy yo*mLAao~ «r- mao 4 O' aoia.-* y mo' y mcom'ymm oiAaoiAaor-Om 
O 4U'm''omcomaomLuma'4oy* , '*r-mo'iA-gr-mo' lAmmiA-gao 4 -gcDin-Ha>irMNjor“vj-Aja'r , **vr cgcM , «-iAivof , -«^uaotnmua)vO'T moo'coaoo'O' 
<X U'coaor-r— yycALA44 mmmmm*-«~*ooo>J'0' -omcor—r-r- y o y a a a 4 4 > 4 - - 4 - mmmmmmm^-g— «^-goooocr*0'0'0' x>auaoaoaor-r— r-r—r— 
JT mmm mmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmrgmm^-g^'-^^ — 



< 

I oO'T>r~yA4m m— «o 0'0»—yA4mm'-«oO' r nr- oiA4mm^oo'''or-yLA4mm^oo''X>r— yiA4mm-*oo'.of s -yiA4 om^oo'aor- OA4mm 
a. oo'0'0'0'0'0'0'0'0'0'aoao''OjQaocoaoaoco<yr-r'-r , - r -r-r--r--r‘'-r— r—yyyy y yy y -o vOcAiAiAininininincr\»A 4 - 4 " nT nT - 4 1 rvj- sj-mmmmmmmm 



<i Huuoouououutuuouuuouuuo'jouuougouuoouuuauouuuouuuoooouoooouououooouauo 

o -* 

u c 

A. oc 



~j (MCNJO^OOrgaOvT OaOAj 40-0 oo^JOgcD<M J^O4C0'0vr4cg r^^-OO O-O'T D'rrg>r'D4vr0 0O>J‘O0'4 J0OOIMMC0 4 >Tvrsrrvlco>rOrgcoO'4-a3 

«/) mm>rao'giA<4-srrgvj-r- ■o^'^oy'^o^^uou^corgm gjmrgr^- g-r-f-o NT'^r^rg D>urgo>0'rcorgcj<-?>ONr coniA'OcM'T - 4 -rgmo o -j- m aj s^mj >o 

> ^»vr gj ^jma'omu''0>'sJOr'-tNj>rO f'i'W -j- f>>Mjau^o -o nOO jj-t rgooo > 4 - ajr'-rgg^u^of'^^'T os^-^ 

uj (Sj'OOrgO'/'O^J'^fNJDO'ocT'rgiA <j ma'too 400 Ajinir>*-< nuo 4 ’^J'' 0 ' 4 ’mof^-' % 4 | AmcgLnvr j'OiAir>oa. , ^NNOr^oo' 4 , m 

O OOO 4<Cinj''-T- o r^-O rg nogpg o 44-^4r''un0^ r ^^rO^GO jMgNmfMM4-OOmO'D^O-n , Pm^J'^^mNcgcMvOC Tj'g^>-*r-^«rin 

■yr (\iooO- <'\jmr^O'rc7'mo>'4’Or s - 4 h jiMn4mM rg.-g ,-\imm <r -nr^ 7)Orjmr-om OCT'M'CO mr- rgoo nc-fa'iAO'^'-"Ofsn3 4'JNmo omor*“ 
<NJO CO O 4rgOCON , ^(r)^Ov’'J00 ^4N-iOO' B^ 0'n4" r HN-<0 coot^ •£ U^>A m rg<>j-g. -tO(J' JOr- s - OcniAinsT 4^mr\j^^^^HQOO> 
(j\Q* aoaoaoaoao^r-r^r^r^p^ -o-o -o ■OvO-0'0'Oinmcnmmmmirunin>r>4 , '4-'f-T'4“'4 - '4‘'r>4->r'r 'T'4'mfAmmmmmmmmmmmmmmmmmmmmmog 



> 

uj cocoh-iriorg^-g'ACT'— ‘msf^T '4“mmrgmmr~oo<?' O' O'O' qjhn o vom<\iC7' -O *0 'jrgor^fA^ >4“ O' A0 0 Jmmmogof s -J' o^ACT'^O'^ >4 - < > 4a' >o<nj iT' r— m 
o miAN>^ 4 , ^omoom>oy'Nio< 0 ''J 00 4 ' ] Of \'00 r aomaO'TO -omo' ■omoN 40 Nmoomo^ ^-gm^rgao'AfNjm^— >r— < io- 4 h nj-4^ 'U 3 > 4 ‘ 
< vf-m '•g— « — «o CT' T' cnh- <; m 'T > 4 -mogrg— •-• oO'^oocd^' sOOinmcrv st g*mmmm*gr^(\j— <ooo U'a'omaj'Dr-^ gg * 0 vOiaialAsT vr nt mmmrgrj 

^ O vO -O-O-O OO'O'OiKK) .0'OOiAiAiALAiAiALAiAiA(AUMAiAiAiAiAiALAi/VALAiAiAiAiA>4- Nr-4'4-4-Nr-rvrs4-'4-Nrg-srg-g--4- s r-4--4--4'>4-Nr-4' 



X Ov7'JOf^-'OiAs4“mrvj'-«0 , >tor^-'OiA^ rr ‘Aj— <OT'00 ,s 'sO^A'J‘fACg-«0«7'vX5f s -’OcA'4“m'\j-HOC'cOf^''OiA>4 , mrg-^0<7'aoh- gjiA •OCT'-Dr^-OcOsTmog 

O. 0 <J'OC 7 'C'^'C f '< 7 ' 7 'C 7 'C' r O r X> t 5 r X)COao' 0 ' 0 ®®NN r, -l s -r s 'r^NN»^N '0 -o ^ O'O'OO >0 'O'OiALAtAiAm.A uAiAiAiA -4“ - 4 “ -4 4444 ' 4 , 4 4 i')mm(<immmm 



<x hOOOOOOOOOOOOOOOO ^>0000000000000000 OOOOOOOOOOOOOO oooooo 000000 00 00 00 000 

o o 

00 • 

a. .x 

* **>**«■**«-**•»»••»•«*•«•«•* o***# 



gjgj r x>ajrg vT g>‘ggjrgrgo^jogoo oO'f OJ-T >rajrg -ro>r rgrg (j'NjOf\J«4 gcuco OO BO 'T '-u O a ■0>0rgoo<rXtcu<rg^;X)OgjJO’U g 

^myy :P'MU'j'N, , ’U 4 g( / J 1 ^^ 004 - 0 ^ oiAOTU'or- 4 '^tr>mr~-^jgjz>s l jg^m'-Jrg g.r > j- N r— <'NJ'-^^ig> 00 >NJ'-D ggjg'j'\y-<Nu\(\i^-<u'\ 

3 ^j- ^j-x -go ’ , nn 4 -f^^>' i iN^mi , JP^®D^n^.mor^uNrg' 7 ' ' 4 D-H 4 ’mu , \-omsrrgu>rg-nmo 043 Nrm'-«r-' 0 ^mco^'>T»or — O -ou^ 

lA^g O'AsO 7*iAO AO OvAocr*- AO'^'A AjgDigg- v/'O'T OrMyN-^r-.NjiAO-OCO O r- gj m g- X) fM vr o'! -O —< O — 4 0 0 O NjmAjr^r^OOOOOO 'J^m— «m<> 
s *^'tiO< , iT\(jivirs.su Qajo f \jm "4 jNrgrgp- i -^j yr\j-H i - 1 \j *o o 'T TOA'OU' r o mm o ' 1 1 imo <Ni <'\j 1 ^ <u v tA O'f^* — ^ 'O m ^g j a NjuAfg'i-o 

^.KTi DO-OnJ - 4\-0 m^O—* 4 3 JJM'|rgn|^uui AJOC 7 ' 0 ''-«' , ^>.C^ ^ <Da^l^J— <~t NrnuVX)i'JN.,n4' vO g r» v A UJ — • -A O' 4 O 4 4 'J- l 04 

,-gys^og <TvOC 7 '<— •'4 r ^0'4X>cg43O'4-C7''4 O' 4 J' O y <\ f^*mo O CM O' lArg -> gj r«)ON ^4 O' -O 4 «-* C4r^ /\mHj'r*inmHOO) 0iA4'M^l4C0r^ gjiAsj- .\l 
4 )’A'A 4 f , »fg.grHO( 7 '( 7 >flur-r- £4iA4 4m^'gMfM'-<^000(M3'004)aONr'r-N>0 0'OiAiAiAiA4 4 nT -4 T ,^1 mm mmm rgrgcg^gcgrg 

mr'^mmmmmmn^AjrgAgog^grgoj'gcgrg.NAjrgrgrgogrgrgcg— ^ ^ w-g -h— 



UJ •gcgin4 4rg^Oco' , ^iAO' 40NiA4 t 14-0 J^maanm 40' 4jao^gs04 O O'AJ O oor- O'— rg 4 -cgrg— «O0'0'rg'04^-<4C0m'0r-.0.0iOf-0'0"A0'O.04' 
o t'Ohoj 4 so aoo 4 O' 4 0*4 o o^r-mo*iA~* r- 4 or- 4 ~.au >04 iAgjcuo.vjiAr»-om>oom«ooiAO k A O' 40 iA-.r-mU ' 43 lM<U 4 0 vO~. 44 -.gjo.A:om 
<i mogomo 4 rg~* 0 'r- 04 mAjocrr- 04 mrgOO' :o oiA 4 NHOO'mr«.p*<nA 4 4 fArgrg^ooO' 0 ' r or-r-r-ooiAiA 4 4 4 r r »mm<NJog<-»~«ooa' aoco 

y 444 mmmmmrgAjrgAjAjrgrg^~«~*^~<~*~«ooooooooo'0'0‘0'crO'0'0'0'0'0'0'0'U'CDa)co«Joac>QDaocox>mtoajcooua3aoaoao3>x>cou>r-r-r- 



X 00 'aor-oiA 4 mAj^oo'aor-o iA 4 mrg ~«o o > aor- <oia 4 mr\j^OO'U 3 r-oiA 4 m<M^oo'aor-oiA 4 mcg~*oO' c *Jr-vOiA 4 mrg^goo'Cur-o^A 4 mrg 
1 1 04 * 4* O' O' O' O' O' O' O' O' au ao qj ao au co co ao au jo r— r- r— r—r—r-r-r— r— r-gj y gj gj -o-o ~o -o 4 J >o umajHiA g> laia lalau ^ 4 4 44444444 m mmm mm mm 

<i -gOO 0000 0000 00 000000 00 00 000 00 0000 00 000 000 000 000 00000 00 000000000 000 0000 

lj m 
u < 



107 



COMPUTER OUTPUT FOUR-THREE 
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OEVIATICNS FOR ALL 
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2000063024. 
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1905584384. 
1860698624. 
1817293056. 
1775303424. 
1734674176. 
1695347968. 
1657269760. 
1620392443. 
1584663552. 
1550044928. 
1516489472. 
1483955968. 
1452407808. 
14218C4544. 
1392116736. 
1363305216. 
1335342592. 
1308194816. 
1281838336. 
1256241920. 
1231379712. 
12C7229696. 
1183765248. 
1160964096. 
1138805504. 
1117269504. 
1096335872. 
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1056201984. 
1036966912. 
1018264832. 
1000 C7 8592 . 
982395648. 
965200384. 
948480000. 
532221440. 
916412928. 
901042176. 
886095104. 
871566592. 
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843713024. 
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817404416. 
804 8066 6 6 • 
792569088. 
780682496. 
769141248. 
757939456. 
747066368. 
736515163. 
72623 8856. 
716374784. 
706766336. 
657462016. 
688456 4 4d. 
679745536. 
671324928. 
663196672. 
655353856. 
647796736. 
640525312. 
633539584. 
626338784. 
620430080. 
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COMPUTER OUTPUT FIVE-TWO 
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PROCEDURE FOR DETERMINING MEANING OF D-STAT FOR SAMPLE WITH 52 OBSERVATIONS 
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COMPUTER PROGRAM TWO 
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START THE FORECAST GENERATION PROCEDURE 
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ROUTINE FOR PLOTTING ALPHA LEVELS VS PODS 

SET NX = NALFA OR NGAMMA DEPENDING OF FORECASTING METHOD 

NX = NAL F A 
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